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ABSTRACT 

The  paper  analyzes  the  h-p  version  of  the  finite  element  method  in 
two  dimensions.  It  proves  that  the  error  |[e||  of  the  method  measured  in 

E  3 

the  energy  norm  decreases  exponentially:  !le[!g  5  Ce-*5  ^  where  N  is  the 

number  of  degrees  of  freedom.  The  exact  solution,  which  is  approximated  by 

2  2 

the  finite  element  method,  is  assumed  to  belong  to  the  space  B  ’  .  This 

p 

space  contains  the  solutions  of  the  problems  of  elliptic  partial  differential 
equations  with  piecewise  analytic  data,  such  as,  when  the  domain  has  corners, 
the  boundary  conditions  and  the  coefficients  of  the  equations  are  piecewise 

analytic,  etc.  Extensive  computational  analysis  with  the  code  PROBE,  shows 
the  practical  effectivity  of  the  h-p  version, 'and  the  applicability 
of  the  theoretical  asymptotic  error  estimates  in  the  range  of  engineering 
computations  and  accuracy. 
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1.  INTRODUCTION 


The  h,  p  and  h-p  versions  are  the  three  basic  approaches  of  the  finite 
element  method  when  the  finite  dimensional  subspace  is  composed  of  piecewise 
polynomials  defined  on  a  partition  of  a  given  domain.  The  convergence  is 
achieved  by  increasing  the  dimension  of  the  finite  element  subspace  in  three 
different  ways. 

The  first  one  is  called  the  h-version  when  the  degree  p  of  the 
polynomials  is  fixed  (at  a  low  value)  and  the  mesh  si2e  h  is  reduced  to 
obtain  the  desired  accuracy.  The  h-version  has  been  investigated  theoretically 
and  practically  for  many  years.  There  are  many  computer  programs  of  the 
h-version  which  are  used  in  engineering.  One  of  the  popular  codes  is  for  example 
NASTRAN  in  its  various  versions,  e.g.  MCS/NASTRAN  and  others  are  ADINA,  ANSYS, 
etc. 

In  the  p-version  the  mesh  is  fixed  and  the  degree  p  of  the  polynomials 
is  increased  in  order  to  reduce  the  approximation  error.  The  development  of 
the  p-version  is  very  recent.  PROBE  is  currently  the  only  commerical  code 
(released  in  1985)  using  the  p  and  h-p  versions. 

The  h-p  version  combines  the  h  and  p-versions. 

The  p-version  was  first  theoretically  studied  in  [  7  ].  The  h-p  version 
was  addressed  in  f3|.  It  was  conjectured  there  that  it  is  possible  to  achieve 
exponential  rate  convergence  in  the  cases  of  practical  importance.  For 
additional  features  of  the  p-version  we  refer  also  to  [  <j].  The  p-version 
in  three  dimensions  was  analyzed  in  I12J,  _£  13 ]  and  for  the  detailed  analysis 
of  the  p  and  h-p  version  in  one  dimension  see  [16  ].  For  the  engineerings 
and  implementational  aspect  of  the  p  and  h-p  version  we  refer  to  [4] . 


It  has  been  proven  that  the  rate  of  convergence  for  the  p-version 
cannot  be  worse  than  that  of  the  h-version  with  a  quasiuniform  mesh  [7]. 

If  the  singularity  of  the  solution  is  located  at  vertices,  then  the  rate  of 
convergence  of  the  p-version  is  at  least  twice  that  of  the  h-version  with 
quasiuniform  mesh.  It  will  he  shown  in  this  paper  that  under  proper  assumptions 
satisfied  usually  in  practice,  the  h-p  version  has  exponential  rate 
of  convergence  with  respect  to  the  number  of  degrees  of  freedom,  while  the 
h  and  p-versions  have  only  a  polynomial  rate  [2]. 

The  singular  behaviour  of  the  solution  of  partial  differential  equations 
of  elliptic  type  is  typically  caused  by  piecewise  smoothness  of  the  input 
data,  by  the  corners  and  edges  of  the  domain,  etc.  Usually,  in  practice  the 
data  are  piecewise  analytic  functions. 

k  £ 

In  Chapter  2  we  introduce  the  spaces  H^’  (ft)  (k2£>0)  which  are  a 
generalization  of  the  weighted  Sobolev  spaces  used  in  [5].  The  main  tool 
of  the  analysis  in  the  present  paper  is  the  countable  normed  space 

0  V  J, 

B  (ft)  which  consists  of  all  functions  u  belonging  to  H0’  (ft)  for  all  k 
p  ,d  P 

and  Hull  .  p  2  Cd^\k-£)!  (£=0,1,2).  It  can  be  shown  £17}  that  the 
Hg’^ft) 


solutions  of  partial  differential  equations  of  elliptic  type  with  piecewise 
analytic  data  belong  to  this  space. 


2-,  2, 


Some  imbedding  inequalities  related  to  the  spaces  H  ’  (ft)  are  derived  in 


Chapter  3. 

The  accuracy  of  the  finite  element  method  reduces  to  an  approximation  problem 


when  the  coercivity  or  the  "inf-sup"  condition  is  satisfied  1 2  j  Ill)  •  We  will 


1  2 

study  the  approximation  (in  the  space  H  )  of  functions  u  f  B0  .(ft)  by 

P ,  a 


the  h-p  version  and  will  show  that  exponential  rate  of  convergence  with 
respect  to  the  number  of  degrees  of  freedom  can  be  achieved. 

In  Chapter  4  we  analyze  approximation,  properties  of  the  polynomials  on  a 


single  square  and  a  parallelogram. 


3 


> 

V 

In  Chapter  5  a  geometric  mesh  on  a  square  and  a  parallegram  is  introduced 
and  the  exponential  rate  of  the  convergence  of  the  h-p  version  is  proven. 

Chapter  6  generalizes  the  results  of  chapter  5  by  introducing  general 
geometric  meshes  composed  of  curvilinear  quadrilaterals  and  triangles. 

The  last  chapter  addresses  the  numerical  results  and  the  performance  of 
the  h-p  versions  by  the  computational  analyzis  of  an  elasticity  problem  by 
the  code  PROBE. 


4 


2.  NOTATION  AND  PRELIMINARIES 


We  shall  denote  integers  by  i,  j,  k,  £,  m  and  n,  and  by  IRi'  and 


IR  we  shall  denote  the  one  and  two  dimensional  Euclidean  space  with 


x  =  (x^jX^)  or  x  =  (x,y).  If  Q  is  any  one  or  two  dimensional  set,  Q 


denotes  its  closure.  By  ft  we  denote  a  polygonal  domain  in  IR  with  the 


boundary  3ft  =*  T ,  the  vertices  A .  ,  1  <  i  <  M  and  T .  »  1  £  i  <  M  the  open 

1  1  M 


edge  of  3ft  linking  and  A^,(Aq  =  ‘  We  nave  3^  =  U  I\,w1iere  is 

i  =  "1 


1=J 


the  closure  of  I\  The  measure  of  the  interior  angle  of  ft  at  A^  is 


denoted  by  w. . 


By  H^ft)  (resp.  Hm(IR2) )  (m>0)  we  denote  the  Sobolev  space  of 

2. 


functions  on  ft  (resp.  IR  )  with  square  integrable  derivatives  of  order 
<m  (ra>0)  furnished  with  the  norm 


u 


Hm(ft) 


.1 


Dau|| 


05(ct|<m  L*"(ft) 


where  a  *  (o^  ,aj,  a±  >  0,  i  =  1,2,  integers,  |ct|  =  +  a2,  and 


<36 

D  u 


3'a*u 


al  2 
3x.  3x0 

X  4. 


u 


ai  a2 
X1  X2 


0  2 

As  usual  H  (ft)  =  L  (ft) .  Further  we  will  use  the  notation 


Hm(ft)  |af=m  p" (ft) 


l  llD““ll2„ 


ana 


|  Dmu | 2 


I  {  |D 

|a|=m 


a  |  2 
u  . 


By  r^(x)  we  shall  denote  the  Euclidean  distance  between  x  and  the 


vertex  A.  of  ft.  Let  8  *  (8, >8, » • • • > 8M)  be  an  M-tuple  of  real  numbers, 
1  1  L  M 


"  5~*  |T»  V^J ; 


i 

I 

ft 

I 


i 


0  <  6  <  1,  IS  is  M,  For  any  integer  k  we  shall  write 

1  mb. 

B  ±  k  =  (B^±k,B2±k, . . .  .B^k)  •  We  denote  <5^  (x)  =  ff  ri'L(x)  and 
M  $.±k  i=1 

Wx)  "  TTriX  (x)- 

By  m  >  l  >  0  (H^,(\ft)  =  K?(ft))  we  denote  the  completion  of 

B  p  B 

the  set  of  all  infinitely  differentiable  functions  under  the  norm 

,2  „  „2 


m  f  k  o  o 

uli^-i(fi) \LJ «|D u<x>l  w*>dx  W51) 

p 


and 


u 


H®’°(ft) 

P 


H.(ft)  k=0 

P 


l  { 

<=n  i 


(*)|2  ^g+k(x>dx 


a  -  o) 


For  m  =  l  *  0,  H^’^(ft)  =  Analogously  as  before 

P  P 


2  -  f 

U  H®,A(fl)  k-A 


l  f|Dku(x)|2  $2+k_jl(x)dx. 


,mi21 


The  space  H  ’  (ft)  was  introduced  in  [ 3  ]  where  its  various  properties  were 
p 

studied.  Let  us  mention  the  following  lemma  proven  in  [5]  which  will  be  needed 
later . 

2  2  0  — 

Lemma  2.1.  H„’  (ft)  c  C  (ft)  with  continuous  imbeddings  i.e., 

p 


sup|u|  <  C||u||  2  ? 
x*ft  Hg  (ft) 


where  C  depends  on  ft  but  is  independent  of  u. 
Let 


0 


(2.1) 

and 


{ u (x)  I  u  €  H™,5,(Q),  Vm  >  l,  0  <  l  <  2} 


6 


(2.2)  Bg(ft)  =  {u(x)  |  u  €  ¥^(Q),  (|  )  DaU{ 2  $g+k.^dx)1/2  <  Cd(k~Z)  (k-£)  ! 

for  |a|  =  k  =  £,£+1,...,  0  <  £  <  2,  d  >  1,  constants 
C  and  d  independent  of  k} 


be  the  countably  normed  spaces  (see  [15]).  For  £  =  0  we  shall  write 

0  £ 

Bg(ft)  instead  of  Bg(ft).  The  functions  in  Bg  (S2)  are  characterized  by 

different  constants  C  and  d.  If  we  would  like  to  emphasize  the  dependence 

on  the  constant  d  we  will  write  B^  (ft)  =  B.  ,  (ft),  etc. 

D,d  p,d 


s 

The  weighted  Sobolev  spaces  H  *  (ft)  of  a  non-integral  s  are  defined  as 

P 


k  l 

the  interpolation  spaces.  Let  8  denote  the  category  of  all  H  ’  (£2) , 

P 


k  =  £,£+1,...  .  8  is  a  sub-category  of  W,  which  denotes  the  category  of 


k  £  k+1  £ 

all  Banach  spaces,  H  ’  (ft)  (k>£>0)  is  a  normed  space,  and  H  *  (ft) 

p  p 


I.  £  _  k  £  k+1  £ 

H  ’  (ft) .  Hence  A  *  (H  ’  (ft),  H  ’  (ft))  is  compatible  couple  in  8.  By 
p  P  p 


the  application  of  the  K-method,  we  define  the  interpolation  space  (see  [10]): 


,„k , £ ..  „k+l,£/  . .  „k+0 ,£,-.  n  .  ,  . 

(Hg  (ft),  Hg  (ft)>Q  ^  -  Hg  ’  (ft),  0  <  0  <  1. 

I 

It  is  easy  to  verify  that  if  u  €  B  .(ft),  then  for  any  k  ^  £ 

p  ,d 

(2.3)  Bull  _  .  i  cjk+"5''-(k+e)1/2  r(k+i-i+6). 

’*<(!) 


* 


»  -  . <■« 


7 


3.  SOME  IMBEDDING  INEQUALITIES 

\ 

We  shall  generalize  some  imbedding  inequalities  proven  in  [2]  [3]  . 

T  be  a  rectangle  (resp.  triangle)  with  vertices  (x^.y^),  1  5  i  5  4  (resp. 

1  5  i  5  3) ;  (x1,y1)  is  supposed  to  be  the  origin.  Let  h  and  h  be 
the  length  of  the  longest  and  shortest  side  of  T,  be  the  measure 

of  the  angles  of  T,  1  <  i  54  (resp.  3).  We  assume  that  there  are  constants  k 
and  oj q  >  0  such  that  for  all  rectangles  (resp.  triangles)  under  consideration 

h/h  5  k  <  co 

and 

0  <  5  0),  5  i  -  w0  <  ’ 

k  M 

Let  Y  be  the  boundary  or  some  sides  of  T,  H  Of)  =  TT«  <Y.),  k  =  0,1 

i=l  i 

and  M  24  (resp.  M  2  3)  where  Yi  ’  s  are  the  sides  of  T,  and  H  (Y  ^ )  is  the 

Sobolev  space  on  Y  . 

i 

2  1 
Lemma  3.1.  Let  u  f  H  (T) ,  and  u  vanishes  at  vertices  of  T,  then  u  f  H  (y) 

and 

(3.1)  iluil 2  <  Chiu!2,  . 

H  (Y)  HZ(T) 

Proof .  Let  S  be  the  standard  square  (0,1)  x  (0,1).  We  obtain  by  the 
standard  imbedding  theorem  (with  C  >  0  independent  of  u) : 


Since  u  vanishes  at  vertices  of  S,  there  is  some  >  0  independent 
of  u  (see  [1]  ,[2]  ,"[!!]) such  that 


and 


Hence 


IU!|  0  5  c ,  M  2 

L"(S)  1  H  (S) 


IU1  1  s  C  |u|  2 

HX(S)  1  H  (S) 


lull2  5  C,  |  u|  2 

rcy)  ‘  h2(s) 


The  usual  scaling  argument  yields  (3.1).  The  proof  for  the  triangle  is 
the  same. 

2  2  B 

Lemma  3.2.  Let  u  $  H„  ’  (T)  with  <*>„=r,0<S<l  and  u  vanishes  at 

B  B 

the  vertices  of  T.  L£  Y  is  a  side  of  T  seoarated  from  origin  then 


□ 


(3.2) 


i  „2.  , 1-2S„  „2,  , 

U  H1 (Y ) ~  h  U  H“ ’ 2 (T) 


Proof:  Let  S  be  the  standard  square  (0,1)  x  (0,1)  and  =  (0,1)  x  (1/2,1) 

Suppose  y  is  the  top  edge  of  S.  We  have  by  applying  Lemma  3.1  (see  also 


// 


Lemma  2.1): 


|u||  i  5  C  |u|  , 

H  O')  H“(S1) 


<  c  u! 


2  2 

Hg  (.S) , 


The  scaling  argument  yields  (3.2) 


Lemma  3.3.  Let  be  the  sides  of  T,  1  S  i  2  4  (resp.  1  <  i  S  3, 

the  right  angle  triangle).  Assume  that  Y^  lies  on  x-axis  and  that  v(x) 
is  a  polynomial  degree  p  on  and  that  v  vanishes  at  the  endpoints 

of  y  .  Then  there  exists  a  polynomial  V(x,y)  of  degree  p  in  x  and 
degree  1  in  y  such  that  V  =  0  on  Y^,  (i  t  1)  and  V  =  v  on  Y^, 

(3.3)  ||V||2  <  Ch||v||2 

H1  (T)  h\'y) 


with  C  independent  of  h,  v  and  T. 

Proof.  Let  T  =  10,1]  x  {0,1]  and  y  =  t(x»0)  I  0  <  x  <  1  }.  Set 

V (x,y)  =  v(x) (1-y). 

Obviously  V(x,l)  =  V(0,v)  =  V(l,y)  =  0  and  V(x,0)  =  v(x).  V(x,y)  is 
a  polynomial  of  degree  p  in  x  and  degree  1  in  y,  and 


Since  v(0)  *  y(l)  =  0,  we  have 


i!v|i22  ^  c|v|2 

LZ(Yl)  Hi(Y1) 

Using  the  scaling  argument  we  obtain  (3.3)  for  the  rectangle. 

The  proof  for  triangle  T  can  be  found  in  [5],  [7]. 

2  2  S 

lemma  3.4.  The  space  H„  (T) ,  <t>„  *  r  ,  0  <  8  <  1  is  compactly  imbedded 
-  P  p 

in  the  space  H^(T). 

2  2 

Proof.  Let  u  £  Hn*  (T)  and  assume  first  that  u(0)  =  0.  Let  (r,0  )  be 

— -  ft 

g 

the  polar  coordinates  and  v  =  rpu.  Then 


6  *  a  S”1 

c  u  a,  +  Br  u  a  ,• 
rQ  e 


o  <  a2  5  2  > 
0  5  a2  5  1  , 


v  2  -  r6u  +  26rB_1ur  +  8(3-l)r6'2u. 

r  r 


Since  u(0)  ■  0  we  have  by  Lemma  8  of  [3  3  (or  Lemma  4.3  of  I  5  ]) 


where  we  denoted 


,  <7  1/2 

B  u  =  (  £  | B^uj  ) 

M=  1 


2.  r 

This  implies  that  v  f  H  (T) .  Suppose  now  that  iu.j  is  a  bounded 

2  ,i=i 

2  2  ft 

sequence  in  H  *  (X)  and  u,(0)  =  0.  Then  v.  =  r  u.  j  =  1,  2 .  is 

p  J  1  J 

2  2 

uniformly  bounded  in  H  (T) .  Because  the  space  H  (T)  is  compactly 
imbedded  in  the  Sobolev  space  W^*q(I)  for  any  1  <  q  <  °°,  ther  exists 


subsequence  denoted  again  by  (v^h^.  which  converges  to  v  f 
w^’q(T)  andv(0)  =  0.  Let  u  *  r  ^  v.  Then 


xn 


~  -8~  p  -8-i~  x 

u=rv-Pr  v  — 
xx  r 


~  -8-  -8-l~  y 

uy  =  r  vx  ”  Sr  v  r 


and  since  v(0)  we  have 


Hence 


ir"s’^n  o  5  cnrVvii  0 

H  (T)  H  (T) 


1°  all  o 

H  (T) 


<  C||r"Wfl 


H°(T) 


by  Holder  inequality 


1 

P 


<  C|l^Spil  0/  , 

n(T) 


!<D1v)q||  % 

HU(T) 


S  C  r 


'BP!I  J  H  i 

rr w/«r\  Ux * 


H  (T)  W  ’’  m(T) 


2q< 


where  —  +  —  *  1,  p  >1  and  Bp  <  1  and  hence  u  f  H  (T). 
P 


Similarily  we  have 


llu  -ull  .  5  Cllv  -v!i  ,  ,  -►  0  as  i 

j  H1  (T).  J  W1,2q(X) 


and  hence  (u^ }  is  a  convergent  sequence  in  H^(T). 


If  u.(0)  ?  0  then  we  define  u.  =  u.—  u.(0).  Then  using  Lemma  2.1 
I  J  J  3 


we 


have 


lu.ll  22  2  C||u  j|  „  -  and  | u  (0)  |  5  C||u  ||  ,  ? 


B 


(T) 


Using  the  first  part  of  the  proof  there  is  a  subsequence,  once  more  denoted 


by  u^  such  that  u^  (0)  A  and  u_j  *»■  u  in  li(T).  This  completes  the  proof. 

n 


Lemma  3.5.  Let  T  be  a  triangle  with  vertices  (x^y^  1  <  i  <  3  and 


=  y-i 


If  u  €  H2,2(T)  with  =  r6,  0  s  g  <  1,  (H2,2(T)  *  H2(T)  for  8=0), 


and  v  is  the 
1  5  i  S  3,  then 


linear  function  interpolating  u  at  vertices  (x^.y^), 


(3. A) 


u-vl!  .  5  Ch3-  ^|u| 

Tt“  / m  \ 


H  (T) 


2  2 
!*■  » *- 


Hg  (T) 


and 


(3.5) 


I  u-v| 


2  2 

Hg,Z(T) 


C|u 


2 

Hg’-(T) 


with  C  independent  of  u  and  T. 


Proof.  If  8*0  then  the  lemma  is  standard  £2].,  £11].  Let  S  be  a  standard 


triangle  { (x,y)  |  0  <  y  <  1-x,  0  <  x  <  l)  with  vertices 

1  <  i  <  3,  and  let  U  f  H^’^(S).  with  $  =  r^  and  V  be  the  linear 

P  P 

interpolation  of  U  (  which  exists  because  of  lemma  2.1).  We  first  prove 
that  for  a  constant  C  independent  of  U 


(3.6) 


Ml2,  ,  5  C(|0| 2  +  l  |U(x  ,y  >|2). 

Hf’Z(S)  H^,Z(S)  i=l 

P  P 


Suppose  that  (3.6)  is  false,  then  there  exists  €  H^’  (S) ,  j 
such  that 


(3.7) 


!fui !! 

J  2  2 

hb  (s) 


(3.8) 


J  |  +  l  0  as  j+®. 

J  ?  9  -{-I  J  1  A 


J  22  i=l  ^ 

H^(S)  1  1 

1  2  2 
Since  H  (S)  is  compactly  imbedded  in  H^’  (S),  by  (Lemma  3.4),  there 

exists  a  subsequence  denoted  again  by  {U.}  which  is  convergent  in  H^(S). 

2  2 

(3.8)  shows  that  {U ^ }  is  Cauchy  sequence  in  H^’  (S)  .  By  the  argument  similar 
to  that  in  the  previous  lemma,  we  can  show  that  there  is  a  subsequence  once 

g  QQ  2 

more  denoted  by  that  {rpU^}  ^  is  a  Cauchy  sequence  in  H  (S), 

r*U  -*■  V  in  H^(S),  and  -*■  U  =  r^  V  ?  H^’^(S)  as  j  -*■  «>.  Since 

|U.  |  -y  9  0,  D^U  =  0  for  | ot |  =  2.  Therefore  U  is  a  linear  function 

3  (S) 

p  2  2  0  — 

of  S.  Because  H0’  (S)  c  c  (S)  by  Lemma  2.1  and  U.(x.  ,y.)  -*-0  as  j  -*■ 
p  j  i  l 

1  2  i  5  3,  U  vanishes  at  vertices  (x^,y^).  Hence  the  linear  function 
U  =  0  on  S,  which  is  a  contradiction  to  (3.7). 


Applying  (3.6)  to  (U  -  V)  we  have 


By  the  standard  scaling  and  mapping  argument,  we  get  (3.6)  and  (3.7).  □ 

Lemma  3.6.  Let  T  be  a  rectangle  with  vertices  (x^y^,  1  <  i  <  4  and 
-  yx  -  0.  If  u  *  Hg’2(T)  with  *  r  ,  0<B<1  (H2,2(T)  = 

H  (T)  for  6  =  0)  and  v  is  the  bilinear  function  interpolating  u  at  vertice 
(Xi.y  ),  1  5  i  .<  4  ,  then 

(3.9)  liu-vjl  <  Ch1"6  |u| 

H  (T)  Hg’^(T) 

and 

(3.10)  Ju-v||  .  0  —  C j  u  1  j  9 

Hg  (T)  Hg  *"(t) 

with  C  independent  of  u  and  T. 

Proof .  Let  S  -  (D-,1)  x  (0,1)  with  vertices  (x^y^,  1  i  c  4, 

U  t  h2,2(S)  with  $  =  r  V  be  the  bilinear  function  interpolating 

c  8 

U  at  'x^v^),  1  S  i  <  4.  can  be  ,  oven  analogously  as  for  (3.6) 

that  for  some  constant  C  independent  of  U 
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(3.11) 


Hull 


2  2 

Hg  (S) 


c(|u|22  2  + 

Hg’  (S) 


4 

l  |u(x  ,y  )| 

i=l 


Let  W  he  a  linear  function  interpolating  U  at  3  vertices  of  S 
other  than  the  origin  and  <f>  he  the  bilinear  function  whic-h  is  equal  to  1 
at  the  origin  and  vanishes  at  other  vertices  of  S.  Let  Z  =  U  -  W,  then  V 
=  W  +  Z(0,0)<}>  and 


[|u-vl! 


2  2 

Hg’Z(S) 


§  llu-wl! 


2  2 

Hp  (S) 


Z (0,0)  |  lull 


2  2 

B 


Ho ’ - (S) 


by  the  imbedding  theorem 


<  c|!u-w|| 

Hg,Z(S)‘ 


In  the  same  way  as  in  Lemma  3.5  we  have 


l|u-w!|  5  C  |  U  |  ? 

Hg  vS)  Hg’^(S) 


which  implies  (3.9)  and  (3.10)  by  the  standard  scaling  argument. 


c 
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4.  PROPERTIES  OF  THE  POLYNOMIAL  APPROXIMATION 


In  this  section  we  shall  analyze  the  approximation  properties  of  the 


polynomials  on  a  square,  rectangle  and  parallelogram. 


Lemma  4.1.  Let  S  =  (-1,1)  x  (-1,1)  be  square  domain  and 
00 

u(x)  *  £  c  P. (x  )p,(x  ),  where  p , (x  )  (resp.  p.(x0))  are  Legendre 

i ,  j=0  1  3  *  il  3  L 


polynomials  on  I  (-1,1)  (resp.  I  «  (-1,1)).  Then 
X1  x2 


f  n  J  0  a  1  9  9  -»  9  (i+a  ) !  (j+a  ) ! 

(4.1)  J  lD  u|  (1-x^  (l-x2)  dx  -  l  c±J  21^1*23+1  (i-a,)!(j-a  )! 


provided  that  the  left  or  right  hand  side  is  finite. 


Proof .  Using  the  basic  properties  of  Legendre  polynomials  we  get  for  m  *  1,2 


rl  (a  )  (a  )  „  ( 

I  m  (xm)pi  *  (1’V 

J  ^  j  mi  m  m 


and  the  lemma  easily  follows. 


2  (i+oO ! 

aH’a^nr  for  i-i 


for  a  >  i,j  or  i  # 
n? 


Lemma  4.2.  If  u  €  H.  (S) ,  k  =  max(k^,k2),  k^,k2  i  2,  then  there  exists  a 

c  i  1 

polynomial  <J>(x  ,x  )  =  l  d.  x  x^  such  that  for  0  £  m  <  2 

OZi±k1  1 

02j<k„ 


-  -  *  -  •  .  v»  w  •  -  ■  r.  r.  i'.  »•  A  ••  •*.  .•  »•  • 


~ T^jv^rv::  /. 1  ■; 1  \  «j  *  ■  «,■  »■  «» r-  ^ 
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(4.2) 


m  o  (k.-s  )!  2 

llD  <“-«ii^s>5  c',v^fcs5T  .1 


s  +1+5. 

|3 _ u  n 2 


l-°  ^,4  h°<s> 


(k  -s  )!  2 

+ - 2_2 -  y 

'•  •  ■  ~  ~  '  i  L 


s  +1+5/ 

1 3  _  u  i,2 


(k2+s2+2-2m) !  ‘  "  ,  •2+l"H0(s, 

3x18x2 


where  s^  is  any  integer,  i  J  s ,  S  K,  (i  =  1,2),  and  C  is  independent  of 


k^.  Moreover 


(4.3) 


u  =  $  at  the  vertices  of  S. 


Proof .  Let  p.(x^)  and  p_.(x^)  the  same  as  the  previous  lemma. 
Since  u  €  H^(S)  we  have  the  following  Fourier  expansion: 

OO 

*  l  ci  .P.U  )p  (x-), 

i , j=0  i  J-  3  z 


2  2 
X1X2 


u  2  Uv-» 
X1X2 


u  2(~1*x2^  * 


X1X2 


u  2(x  ,-l) 
X1 


u  2(-l,x2) 
X2 


l  a  p  (x  ) 
i=0  11 


^  b  p  (x  ) , 

j=0  jV  2 


I  Vi^i' 

1=0 


-  I  Yj(*2> 


3=0 


Set 


Xi ^xi»x2> 


I  ci  -jWW 

0<i<k  -2  1,3  1  1  3  2 


0<j-k2-2 


X2(x1> 


0<i<k1-2 


aiPi(xi)’  X2(xx) 


A 

■  J_iX2<5i)dCi  +  uV2<-l.-l). 


X3(x2) 


y  b  p.(x  ),  x_(x„) 

<i<tL-2  3  3  2  32 


02j 


X4(xx) 


l  V^V’  W 

0<i<k  -2  1  4 

1 


X5(x2) 


l  X5(x2)  = 


05j<k2-2 


J_x  X3(?2)d?2  +  Ux1x2("1,"1)’ 
rXi 

J.L  X4(«l)d€l+ttx1<-1-“1>* 

rx2 

j.,  x5a2)d^2  +  UX2<-1*"1)- 


Let 


(j)(x)  =  <f>1  +  4>2, 


where 


x  x  £  £ 

<4'4>  *1  '  J-!  Ll  Li  Li  Xl^l’VVW«2  +  <1+V  {_!  VVS 


r:ii 


+  j ,  h(hHh 

-1 


(A. 5)  <t>2  =  (l+x^l  X3(C2)dC2  + 


rX2 


X,(5  JdC,  +u(-l,-l)  +  (1+x.  )  (l+x,)u  (-1,-1) 

5  2  2  1  2  x^x2 


Hie  degrees  of  4> ^  (i  =  1,2)  in  x^  end  x2  are  at  most  k^  and  k2 


respectively.  We  can  readily  verify  that 
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x  x_ 

f  1  f  2  f  x  f  2 


-  j.x  j.,  j_.  j.L  “A2«i:vwii5!t(iwi)  L-W1,5**5* 

1  2 

x.  x  x 

r  i  f  ^  2 

+  (l+x  )  u  u  (5  ,-l)d?+j  U  (£  ,-l)dC  +  U  (-1,£ 2)d^  +  u(-l,-l) 
*■  J-i  xi  x2  L  1J_i  xi  ■*-  x  _i  x2 


(l+x  )(l+x,,)u  (-1,-1). 

-L  a  X  |  x  2 


U  =  U1  +  U2 


x_  x0  C. 

rlr2flf2 


(4.6)  u.  = 


I  <  I  +  I  He  P  (£  )P  <e  ))d5  d5,d5  d5, 

J-l  J-i  J-i  J-i  0Sl*k  -2  i>k  -1  1  1  1  1  " 

JL  1 


05j5k2-2  j>0 


+  (1+X2}  J.x  V^l’^l  +J-1 


"Ddq, 


x  x  C,  Co 

r  1  r  2  f  jl  r  2 


(4.7)  u£  = 


f  1  f  M  A  f  *  r  ~  -  ~  ~ 

i  I  c  iPi(VP1(Vd^-J^2 

J-l  J-l  J-l  J-l  02i<k  -2  1,3  x  1  J  1  1  • 


1 

jSk9-l 


X2  C  2 

+  (l+x  ) [  u  (-1,?  )dC0  +  [  u  (-1,£  )d£  +  u (—1 »— 1) 

i  J  x]/^2  ^  J  x2  ^  ^ 


-  (l+x  ) (l+x  )u  (-1,-1); 
1  2  x^x 


u  —  <{)  ■*  u  -  <|>  +  u2  -  <j>2 


F1  +  F2  +  V  U2  ’  *2  =  F4  +  F5  +  V 


we  get 


i-0 

X  x  £  B, 

r  1  r  2  r  1  t  2 

I  c  p  (C1)p.(?2)dC1dC  dc  dr  , 

-1  -1  J  -1  -1  0515k  -2  l»J  x  1  J  2  1212 


j*k2-l 


rX2  A 


=  (i+x  )  [  [  y  aipi(^i)d?i’ 

j_i  j_i  i>k  -i  111  x 


x2  B, 2 

d+V  |  f  I  bp  C)d?dr  , 
J-l  J-1  i>k„-l  ]]  ‘  1  Z 


fXl  A 

"  I  diPi(Cl)d^ld?l’ 

J-l  J-l  i>k  -1  111  1  1 

x2  ^2 

-  f  1  l  e  p /C  )d?,dC2. 

J-l  J-l  1>k„-l  332  ‘  2 


It  can  be  seen  that 


x0  Co 

f  2  f  2 


J-l  j-l  Jrl  Cl,3Pi<Xl)PJ<?2)d?2d52 

05j<“ 

1  /  P^i  .o(X9)_P  -  (Xo> 

y  c  p  (x  )  — =— i  (  —  Ji± — 2 — i — - — 

i>k  -1  i*i 1 1  (2j+1)  "  2j+3 

05  j<^ 


Pi(x2)-P1-2(X2)\ 

2J-1  ) 
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,1IV  <  c  y  _ * _ i _ 

'  2  'I  0  ~  C  *  .  21+1  (2j+l)  ......  A 

SXj^  J  (2j+l) 


(k-s)< 


.  — — — — —  y  |c  r  _ £ 

c:  +s  -2)  1  A  1  i-  -i  I  2i4 

11  )m  i>s  -1  l»J  ^ 


2  2  (i+s1-l) ! 

2i+l  2j+l  (i-(s  -1>) ! 


by  Lemma  A.l 


(k-s)!  Sl+3 
1  ||_j£ _ U _ || 

W^'  V\  2  H  (S)  ' 

3xx  3x2 


In  the  same  way  we  have  with  C  independent  of  k 


(4.3)  l|D"Fll^0(s)  -C  (kl+i1+2-2.)!  0S" 


s  +3 

3  1  ii  2 


ox^dx^ 


Next  we  estimate  F2 .  Since 


“  (  I  a.p.Cx^Hl^) 


i>k  -1 


2„2  <  y  2 

2"  0(1  )  ~  C  L  a-« 

3x^  r^x^  i>k1~l 


2  _i_  <  (ki~si)! 
i  2i+l  “  ^  (k.+s  -2) 


(i+s  -1) ! 


1  1  _  y  2  — 2 —  1  _ 

,+s  -2)1  l  .  ai  2i+l  (i-(s  -2))! 
1  1  i>s  -1  1 


<  c  (krsi)!  11 2 

<kl4v2) !  ,  Sl+1a  H°(I  ) 

3x^  8X2  x^ 
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we  obtain 


0  r  o  2 

Hr"  o 

h  (s) 


In  a  similar  wav  we  have 


(k, -s. ) ! 


'Sx  11  0 
9X1  H  (S) 


If2II  o  s  c 

H  (S) 


(k1+s1-2) ! 

(k  -s.)! 

X  X  II. 

■  H2 - - — 

i!  s  +1 

3x11  3x2 

s  +2 

3  u  „2 

(k+s .)!  H 

s  +1  II  ( 

1  1 

3x]_1  3x2  h 

(k  -s)! 

,sl+2 

1  1 

r\~  \  i 

II A  .  “  ! 

H  (S) 


3*,1  3,2H(S) 


32F  8F, 

It  is  obvious  that  — j-  =  0  and  ||— — 1| 

3x2  3x2  H  (S) 

for  0  5  m  2  2 


-  c!|F„|j  o  •  Therefore  for 
2  HJ(S) 


(4.10) 


m  9  (k-  ^  ’  n  -3-  o 

Id  f  II2  c  r  _ i  F _  |i _ 3 _ y _ 11^ 

2  H°(S)'  Vs1+2-2->!  "  V2  H°(S)  * 

3xx  9x2 


s  +2 

3  u _ ii  2 


Analogously  we  have  for  0  s  m  <  2 

9  (k  - 

(4.11)  ||D  F  IP  5  C  TTTTT- 

5  a  (s) 


(VS2)! 


(k1+s2+2-2m)! 


s„+2 

3  u  ,i  2 


S2+1  V(S)  ’ 

3Xl3x22  (S) 


(4.12) 


(Jr  _S  )  t  S.+l 

l|DmF  I!2  <  C  1  1  i|3___Jl|| 2 

3"h0(S)‘  Ck1+.1+2-2-)l  'I  VlV(s)- 

3  x1 


fk  —  s  )  * 

ra  2  so '  •  %  -  ,,  o 

|r\IU-p  ii ^  <  r  _ **  ^  i;<7  u.i*- 

SH°(S)~  (ki+3,-2-21) I  '  s+10 

a  x 


(4.13) 


Hence  ^sing  (4.8),  (4.10)  and  (4.12)  we  get 


(4.14) 


(k  -s  ) '  2  Sj+l+£ 

||Dm(u  -♦  )!|20  <  C  ^.2  ,  l  II3  .  +1  u  II 2U  ,  osi;! 

1  1  H  (S)  U1  S1  Z  1=0  „  1\Z  H°(S) 

3x  9x2 


and  by  (4.9),  (4.11)  and  (4.13) 


(4.15) 


,2  (VS2)! 


2  S2+1+5, 

r  l  iihrns«2n  • 


I|D  (U2'*2>"h°(S)  1  C  (k2+S2+2-2”)!  Jo  "Tv5^ 'As) 

9x1 


(4.14)  and  (4.15)  yield  (4.2).  From  the  orthogonality  of  Legendre  polynomials 
F1(±1,x2)  =  FA(Xl,±l)  =  F2(±l,x2)  =  F5(x1>±1)  =  F  (±1)  -  F6(±l)  = 

which  implies  (4.3).  | 

By  the  scaling  argument  we  obtain  immediately  the  following  lemma 

Lemma  4.3.  Let  ft  =  (a,b)  x  (c,d)  with  h^  =  (b-a)  and  h2  =  (d-c) . 
k+3 

If  u  €  H  (ft),  k  =  max(k^,k2) ,k^,k2  >  2,  then  there  is  a  polynomial 

<{>(x)  =  J  c  .  .x^x^  such  that  for  any  integers  15  s.  -  k.  ,  i  »  1,2,  am 
05i5k1  1,J  1  1  1  x 

0<j5k2 
0  <  m  <  2 


(4.16) 


/,  \  |  t.  2(s  +1  0  Sl+l+£  , 

m  2  -2mf  (kl~91)  *  \  1  f  ..3  '  u  „2  \2 Z 

J  H°ffi>  -  -  <  W2'2”’ '  *  ’  l0  i  s1+l  »H0(n)(  2 

d*2 


(VS2)!  h2  2(S2+1>  2  ^2+1+\  2 

+  (k„+s  +2-2m) !  £  I!  r+T'l 


T  ^  h 

1 - Hr 2  (-V) 

n  S„+1.0  2  > 


where  h  =  minCh^h^)  and  C  is  independent  of  k.  Moreover 


(4.17) 


u  =  cf>  at  the  vertices  of  ft. 


Theorem  4.1.  Let  ft_L  =  (a,b)  x  (c,d)  cc  ft  with  ^  =  (b.-a)  <  XjT  and 
^2  =  (d-c)  <  ^2r0*  ^i  -  (i  =  1.2)  and  r^  be  the  distance  between 

origin  and  ft  .  Assume  that  h/h  =  maxO^.h^/minO^.i^)  <  <.  If 
k+3  2  8 

u  €  H  ’  (ft)  with  $=r,  0  <  g  <  1,  then  there  exists  a  polynomial 


<n  X  i 

4>  =  l  c  X  X7  on  ft.  with  2  <k.,  k  <k  such  that  for  0  <  m  2  2 
0<i<k  1  z  1  1  z 

0<jSc2 


(4.18)  |!Dm(u-<J>) 


_  9  9/9  Q\  ^  r(k.~S.+l)  X.  2S . 

'  (u-<f,)li2„  <  cr2^2“ra“^  Y  - 2 _ 2 _ (-J-)  3  ||u|| 

V  11  n  ~  n  L  P/l,  J-o  4-9_9™N  v  O'  llUll 


1  0 

H (ft^) 


.L.  r(k.+s.+3-2m)v  2  11  11  s.+3,2 

J"1  J  J  u  2  < 


s .+3,2 

where  s  is  any  real  number  1  <  s.  <  k.,  j  =  1,2,  and  H„J  (ft)  is 

j  ~  2  ^  J  p  ' 

k.+2,2  k  +3,2 

the  interpolation  space  (H  J  (ft),  H  2  (ft)).  for  some  integer 

o  p  Wj »°° 

k.  *  s  +  1  -  8,  S  k  ,  0  <  0  ,  <  1  and  C  is  independent  of  k,  but 

J  J  J  J  J 

dependent  on  X  j  and  s^ ,  j  =  1 , 2 .  Moreover 


(4.19) 


u  =  at  the  vertices  of  ft 


Proof .  Applying  Lemma  4.3,  by  (4.14)  and  (4.15)  we  have  for  the  integer  k 


_  2  2™  -k. ) !  h.  2(k  +1)  2  , 

l|Dm(u  -*  )||20  <Ch2”  - i-i— <-i>  n»E+J  II  0  ^ 

1  1  H  (Q  )  (k+k+2-2m)!  2  £=0  V1  l  H(fi  )  2 


(k  -ic  )  *  X  2^i 

-  Cr2(2-m-S)  U1  V*  1  ,<1,2 

cW2-2*),  A  3>  (B) 

P 


and  similarly 


|Dn(u  -*  )||20  <  Cr2(2-n-6)  — (^  k2>  ' -  (~)  2||u||\  +3  2 

H«)  (k2-V2-2„),  «g2  («) 


X  «  =  D“(u  -$  ),  3  -  1,2. 


k.+3,2 


Thus  T,  is  an  operator:  H  2  (Q)  -*•  )  and  the  norm  of  the  operator 

3  8  1 

is  bounded 

fk  -lc  3  •  \  ^4 

III  ll2~  <  Cr2(2-m-8)  ■„:kl  Y'  A  2 

i  k.+3,2  ~  0  „  o  x.  '  l’ 

H3j  (kj  kj+2-2m)  • 

s.+3,2 

If  s.  =k.  -1+0.,O<0.  <1,  H  J  (G)  is  defined  as  interpolation 
3  1  3  3  3 

space  [10] 


s .+3,2 


k.+2,2  k.+3,3 


H„]  (P)  =  (H  ->  <G),  H„J  (fl)). 


s .+3,2 

T.  is  also  linear  and  continuous  operator:  H  2  (£2)  )  and 

3  3  1 


|Tj  11  8  +3,2  "  l|Tj!i  k  +Vt2  ,||Tjllk  +3,2 

h3j  w),Hra(a1)  h  J  (g^h"^)  h3j  (G), 


2(0  m  (k.-k.+l)1-8 

<  Cr2(2-m-S)  f_Jx  _ j _ 3 _ \ _ ] _ 

u  (k.+k.+2-2m) !  (k.+k+2-2m) 

3  3  3  j 


For  any  integer  n  >  0  and  0  <  9  <  1,  there  are  constants  C^,  >  0 

independent  of  n  such  that  [18,  p.937] 


(4 .20 ) 


Cir(n+l+0)  <  n!  (n+1)  <  C^n+l+Q) 


Hence  we  have  for  j  =1,2 


o  90  q\  X.  zsj  I\'k  -s  +1) 

n  i| 2  <  Cr2(2-m-8)  /_!)  - 3_J - 

Lj»  s  +3,2  "  0  L2'  r(k  +s  +3-2fl) 


„  OO  AN  2  Kk.-S.+l)  X.  j  , 

Id  <u-4>) ||  0  <  Cr  l  r(k  +s  +3_2m)(  2_)  ^  s.+3,2 

H  (Q1)  J-l  3  3  H  3  (Q) 

p 


(4.19)  follows  now  from  (4.17). 


Let  ft  be  a  parallelogram  with  interior  angle  u,  0  <  w  <  ^  and  two 
sides  a  and  b  (see  Fig.  4.1).  We  introduce  the  mapping 


Fig.  4.1 


Parallelogram 


x^  =  X1~X2  cot  w 


x2  =  x^/ sin  to 


and  have  analogous  lemma  and  theorem  for  a  parallelogram  domain. 


Lemma  4.4.  Let  ft  be  a  parallelogram  shown  in  Fig.  4.1.  If 
k+3 

u  €  H  (ft),  k  =  max(k^,k2),  k^,k  ^2,  then  there  exists  a  polynomial 

<f>(xi  ,x  )  =  1  c  .xl’Xo  such  that  for  0  5  m  2  2 

0212^  1,2  1  Z 

0<j2k„ 


(4.21)  ||D“(u-*)|r0 

H  (ft) 


(k  -S  )! 

SC(min(a,b) )  *  (kj+Sj+2-2111)  ! 


2(s  +1)  21 

f*  l  l  !lu  s  +l+o  -  II  0  <f> 

2  £=0  l'=0  1  1  l-V  H  (ft)  2 

X1  2 


(k  -s  )!  2(s+l)  2  S2+1  a 

+  (k  +s  +2 -2m)  '  °  11  ilu  s  +1-2,'  II  ^  ) 

U2  2  m) •  £=0  V=0  l+l'  2  0/o. 

xx  x2  H  (ft) 


where  X^  *  ^xi"x2  cos  ^2  =  x2’  si  is  inte§er  -k^  i  =  1,2  and  C  is 


independent  of  k.  Moreover 


(4.22) 


u  *  <j>  at  the  vertices  of  ft  . 


Theorem  4.2.  Let  ft^  be  parallelogram  ccft  with  interior  angle  to, 

IT 

0  <  ui  <  ~ ,  and  let  a  and  b  denote  the  lengths  of  its  two  sides, 

a  2  a°d  b  -  ^2r0  wkere  r^  is  the  distance  between  ft^  and  the 

k+3  2 

origin  where  the  vertex  of  ft  is  located.  If  u  (H.  ’  (ft)  with 

p 

$R  *  r^,  0  2  6  <  1  then  there  exists  a  polynomial  <j>  =  £  c.  .  xjxj, 

p  n<  i  <!/■  » 3  J  i 


2  2  k^,k2  <  k  such  that 


0<i<kL 

0<j<kR 


(4.2  3)  ||Dm(iH>)l|20 


H  (fl) 


o/o  2  r(k  -s  +1)  A,  2s  , 

<  rr2(2-m-8)  c  - j__j - , — 2_.)  J  i|u|| 2 

-  Cr0  r(k.+s  +3-2m)(  2-y  Ul1  s  +3,2 

J  J  *  H  J  I 


where  s.  is  anv  real  number  1  S  s .  S  k ,  (j  =  1 ,2) ,  C  is  some  constant 
2  11  si+3>2 

independent  of  k  but  dependent  of  angle  ,  (Q)  is  the  interpolation 

lc  +2,2  k.+3 ,2 

space  (H  2  (Q) ,  E2  (£2))fl  .  0  i  0  S  1,  for  integer  k.  <  k.  such 

8  p  p,00  1  1 

that  s.  =  k  -1+0.  Moreover 
1  1 


(4.24) 


u  =  <}>  at  the  vertices  of  (2 


5.  APPROXIMATION  PROPERTIES  OF  PIECEWISE  POLYNOMICAL  ON  A  SQUARE 


AND  A  PARALLELOGRAM 


Let  ft  be  a  square  or  a  parallelogram  and  u  (  ^(ft)  with  the  singularites 

in  the  vertices  of  the  domain.  We  shall  show  that  h-p  version  leads  to 
an  exponential  rate  of  convergence  if  the  mesh  and  the  distribution  of 
the  degree  are  properly  chosen. 

2 

Let  us  investigate  the  approximation  of  the  function  u  €  B„  (ft'* 

p,d 

g 

where  ft  *  (0,1)  *  (0,1)  and  *  r.  Let  a  mesh  on  ft  be  as  follows: 
let  Xq  *  =  0,  and  x.  ■  y.  =  an+^  ^ ,  1  ^  j  £  n  +  1  for  0  <  a  <  1  and 


“i.j  ‘  fcj-i’V  *  <3rj-i,yj> 

S2,J  "  X 

n3,3  =  <0’x3-l)  *  Vj.i.yj) 


for  1  £  j  £  n  +  1  , 


for  2  £  j  £  n  +  1 


for  2  £  j  £  n  +  1  . 


The  nodal  points  which  are  marked  *  in  the  Fig.  5.1  are  classified  as  irregular 

nodal  points  and  the  others  as  regular  nodal  points.  The  rectangulation 

{ft  ,  i  =  1  for  j  =  l,  l£i£3  for  1  <_  j  £  n  +  1  }  is  denoted  bv  ft*n. 

J  0 

k 

If  we  divide  ft.  .  for  i  2  2,  j  >  2  into  three  triangles  ft.  ., 
ij  s  i  *  J 

1  £  k  £  3  as  shown  in  Fig.  5.1  we  get  a  mesh  on  £* s  (ft.  , ,  for  i  ■  1, 

1  *  J 

1  £  j  £n  +  1,  2  <  i  £  3,  j  -  2,  ftj^  for  1  £  k  <  3,  2  £  i  £  3,  2  <  j  £  n  +  1}  . 
We  denote  this  mesh  by  ft^  and  call  it  a  geometric  mesh. 


m 


^,3  |^2,4 


Figure  5.1 

Geometric  Mesh  on  a  Square 


Let  S?’2’V)  c  H1^)  be  subspace  of  functions  which  restricted  to 

a 

ak  ,,  lSk<3,  2<i  <3,  2<j<  a  +  1  (resp.  Q  i  »  1,  lSjSn  +  l, 
3-*  J  1  »  J 

l<.i<  3,  j  =  2)  are  polynomials  of  degree  <  p.  .  (resp.  p.  )  in  x  and 

1  f  J  1  9  j 

of  degree  -  q.  .  (resp.  q.  .)  in  y.  By  P  and  (£  we  denote  the  degree 
i*  J  1 » J  “ 

vectors  P  -  (p^.)  =  <Pia.P1>2.P2,2»  *  * '  ’P3,n+1)  and  (qi,j}  = 

(ql,l,ql,2’q2,2’*"’q3,n+l)* 


2  8 

Lemma  5.1.  Assume  that  H  =  (0,1)  *  (1,0),  u  f  BR  ^(Q)  and  *g  =  r  .  Then 

there  exists  <J>(x,y)  f  (fl”)  ,  0  <  o  <  1,  p^  ^  s  q^  ^  =  1,  pk  .  -  2, 

q’f  .  =  q.  .  £  2,  1  <  k  <  3,  j  Z  2,  and  p,  . ,  q.  .  are  non-decreasing  with  j 
i,j  x,j  i>3 


such  that 


l<i<3,2<j<n+l  j_1  (pi,j+Si.j“2) !  '  2  ' 


V  ‘“i.j’ 


2(l-6)  (ql,3  ‘l.J1-  p^L  1,J|,  „2  .  } 

!j-l  (q  +t,  -2)!  \  2  I  Wl  t.  ,+3,2  >  } 

,J  ’2  to,  .) 

P  1  *j 


where  <l>-|  -i  “  <J>(0  is  bilinear  interpolation  of  u  at  vertices  of 

at-i  ■* 


n 

“i  i  » 


s.  .  and  t.  .  are  real  numbers,  1  <  s.  .  5  p.  .,1st.  .5q.  . ,  X,  .= 
i.J  1.3  i»J  i .3  i.j  i+l  1.3 

^2.  j  =  (1*  o)/(/2  o),  —  v3  j  —  1*2  j  =  X^  j  =  1  for  2  £  i  S  n  +  1*  and 


(pi,j  "  si,j)!  *  r<pi,j  '  *i,j  +  1,1  etc- 


Proof.  Applying  Theorem  4.1  on  each  Q.  . ,  1  <  i  <  3,  2  <  j  <  n  +  1,  there 

i » J 

exists  a  polynomial  <f> .  .  with  degree  p.  .  in  x  and  q,  ,  in  y  such  that 

i»J  i»3  i.j 

for  0  s  m  2  2 


(5.2) 


!|Dm(u-<J,  ,)||2q 

’3  «  (%.j) 


<  C  x2(?_m_B) 


2s 

(Pj  , )  •  /  X.  •  .  v  i.j  2 

p/r _ 1.1  l»j _  1  .J,  Hull 

j-1  l(p.  ,+s  +2-2m) !  \  2  /  ,|U"  s  +3,2 

’2  H„  to,  .) 


(qi>^~<:i.j) !  /Vi.j 


(qi,j+ti.j+2_2m) !  ’  2 


t  ,+3,2 

H6  3 
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On  the  regular  nodal  point  of  ft*n  the  polynomial  <j> .  . (x,y)  coincides 

O'.  1  ♦  J 

with  u(x,y),  We  need  to  adjust  <j> .  .  on  ft.  .  for  the  coincidence 

on  the  irregular  nodal  point. 

We  show  the  treatment  on  ft,,  .  (j  >2).  We  denote  <j>.  . (x,y)  on 

k  k  k 

ft~  .  by  <j>  (x,y),  1  5  k  5  3,  and  let  v.  .  be  linear  function  on 

^»3 


ft-  .  such  that 
2,3 


k  k 

v  .  (x._,  ,y._2)  =  u^._1,y._2)  -  <t>2  a  » 

^3  3  3  3  J  » J  3  3 


k  k 

v_  .  =  0  on  the  other  two  vertices  of  ft-  . . 
2,3  2,3 


Then  setting 


-rk  ,k  .  k 

?2,J  -<-2,3  +V2,j> 


on  ft,  . ,  IS  k  53 
~  *  J 


we  have 


k  k 

u  =  3L  .  on  all  vertices  of  ft„  . , 

y2,3  2 , j 


(xj-i»yj_,) 


(xj-pV2  )' 


Uj.y,.,) 


fl2,i 


ft2,i 


(xj_i,0) 


(*i,0) 


1  <  k  <  3. 


Figure  5.2 


Partition  of  a  Nondiagonal  Element 


Applying  Lemma  3.5  we  have  for  i  2  2,  j  >2 


a 


(5.3) 


(5.4) 


!|u^J22  k  ^  l|u-»J  -v\  ||22  k  iClM.J  I2, 

'3  HZ(fi^  )  1,J  ’J  tT(fr  .)  i,j  H  (ft*  .), 

1  *  J  ^ » J  1 » 3 


.  k  .,2 


,k  ,  2 


Mi  jli  !  k  5  Ch  |u-4»  |  2  k 

■L’J  h-vo*  \ 


ror  .) 

i.J 


IT  (ft?  .) 
J 


fk 

Functions  (J)  .  (i  =  1,  or  j  <  2)  and  $ .  .  (i  >  2,  j  >  2)  should  be 

J  J 

further  adjusted  for  the  continuity  on  the  common  edge  of  two  triangles  or 
rectangles . 

There  are  six  basic  cases  for  i  ■  2  shown  in  Fig.  5.3,  and  another 
similar  six  cases  for  i  =»  3,  which  we  have  to  treat. 
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Since  $9.-$9.=v9.-v  is  a  linear  function,  w  vanishes  on  Y . 
z»3  z»3  z»3  z>3 

1  3 

Hence  the  piecewise  polynomial  is  already  continuous  on  S3.  .  U  fl9 

z » J  z  ,3 

2 

(j  >  2).  In  Case  2,  the  functions  <j>  .  ,  and  $  .  coincide  with  u  at 

1 » 1  1  3 , j 

2  „2 

the  vertices  of  .  -  and  ft-  . ,  w  =  ($..  .  . )l  vanishes  at  the 

-*- » 3  z»3  i»3_i  z»3  J 

2 

endpoints  of  common  side  y  of  ft  .  1  and  ft-  . .  By  assumption  w  is  a 

1 » 3 “1  z » 3 

polynomial  of  degree  5q  .  in  y.  By  Lemma  3.3  there  exists  a  polynomial 

z » 3 

A 

w  of  degree  q9  .  in  y  and  degree  1  in  x  such  that  for  0  5  y  5  y.  .  , 
z»3  3“L 

~  2 
w(x.  9,y)  =  w(y),  it  vanishes  on  other  sides  of  ft9  and 
3-z  z  >  3 


(5.5) 


Iwli  .  ,  —  Ch„  .llwll" 

H1®,  )  2,3  H1  (y ) 

L  9  J 


5  Ch  .{llj2  .  -  ul|2  +  l!d>  .  -  ull2,  ) 

2.31  2.3  Hl(y)  1.3-1  Hl(y)l 


by  Lemma  3.1 


<  C{h 


2  ,,rr2 


2  1  ^2  i  '  U"  2  2  +  h2  i  h2  1-1  ii<J)l  1-1  ‘  u!l  2 

2,3  2,3  H2(ft2  .)  2,3  20  1  1,3  1  H2 

^  >  J 


(Qi,3-i} 


<  c{h2  ll$2  i  -  ui|2  2  +  h2  j-l^  1  1-1  ‘  u!|22 

2,3  2,3  H2(ft?  .)  2,3  1  1,3  1  H2(0,  ,  . 

z»3  —  »3 


Setting 


2  -2  9 
>P9  .  =  ,  -  w  in  S2_  .  , 


^1  .  =  t{)  .  in  ft.  .  . 

l.J-1  l,j-l  1,3-1 


we  have 


(5.6)  Hu  -^J2!  2  -  2{\\u-Z  2  +M\  2  } 

2-  ulrnz  ^  z>3  H  (ft  t  .)  H  (ft,  .) 

*•  ’  3  2  *  2 


-2  ||  2 


H  (ft,  .) 

^  >  J 


-  {!!u_fP?  J21  2  +  h2  il!u"^2  1^2  +hi  i  ill  u~\  i-1^2  } 

“’3  H^flf  .)  Z'2  l'2  H2(ft,  ,)  1»J”1  1,3  1  HZ(^  ,  ) 

2-  *  J  +-  9  J  ■L^J  -L 


r2  ||2 


'i  2 


by  (5.3)  and  (5.4) 


-  Cth2,J  l“-02,-'2 


•i  I  2  k  +  ^1  i-l^U  ^1  i-ll  2  ^ 

3  HZ(flT  .)  ,J  1,3  llHZ(fi  ) 

2>J  -L'J~X 


and 


(5.7) 


lu-*i  j-i11  i  =  !!  u~$i  i 

,J  1  ,  .)  n  ( ft.  ,  .) 

1  1  J  “'  X  1  *  J  X 


In  Case  6  (j  =  2)  ^  ^  is  the  bilinear  interpolation  of  u  at 

vertices  of  ft  The  common  edge  of  ft..  .  and  ft  is  separated  from 

-1*1  -L  ,  1  ^  y  ^ 


the  origin  and  w  =  (4>2  2~^l  vanishes  at  endpoint  of  y.  By  Lemma  3.3 

there  exists  polynomial  of  degree  2  in  y  and  degree  1  in  x  such  that 

A  A 

w|  =  w  and  w  vanishes  on  the  other  edge  of  ft,  9 


and 


IlSil  ,  <  Ch  Jlwll" 

H  <!22,2>  ’  H  (V> 


by  Lemma  3 . 1  and  Lemma  3 . 2 


<  Ch  -(h  ,||w-4>  ||  _ 

l, l  l, l 


+  <‘?ell“-+,  /2  2 

2>2)  1,1  ’  BJ’  (»! 


Letting 


we  have 


2  2 

-  C(h-  J|w-<|>-  _||  -  >  1*1  1  1 1I  2  9 

2’2  2’2  H2(n2>2)  1}1  1’1  H2’2(^1>: 


r 


lb  m  4>  -  W 

r2,2  2,2 


«  ib  =  <t> 

^  1,1  1,1 


in  ^ 


2  2 


in  Qia, 


(5.8) 


u-ip.  A2,  <  !|u-<f2  /,  +  llwll2! 

’  »  <122,2)  ’  H  <H2,2)  H  (SJ2,2) 


<  c{»u-*  a2  +  h2  >-*,  ,r2  +  hrr°/||u'h,i 

’  H  (!I2,2)  ’  ’  H2(!J,  ,) 


2(1-6) 


(^2,2^ 


and  by  Lemma  3.6 


lU-'A 


s  Ch2(}-®»u-. 


(5.9) 


9  9 


After  treating  all  the  cases  we  obtain  modified  V  ~  ,  i  >  1,  j  >  2  and 

J 


modified  ^ .  (i  =  1,2).  Define 


.  on  .  for  1<  k<  3,  2<  i<  3,  3<  j5  n+1, 
i>J  i,J  J  ’ 


.  on  ft  .  for  i  =  1,  1  5  j  5  n+2  and  1  <  i  <  3,  j  ■  2. 

*  J-5  J  1  J  J 


Hence  $  (  and 


(5.10) 


u-C,  ^  C{hJ(J"6)||u-4)  || 

H±(Q)  1’1  HX^Cn,  .) 

p  1,1 


+  l  Clf  u-4>  II2.  +  h2  i|u-4>  ||20  )}. 

l<i<3,l<j<n+l  H  (ft.  .)  1,J  1,J  H  (ft,  .) 


(3.5),  (5.2)  and  (5.10)  yield  (5.1) 


Theorem  5.1.  Let  ft  =  (0,1)  x  (0,1)  and  ft^  be  the  geometric  mesh  on  ft 

2  8 
shown  in  Fig.  5.1.  If  u  t  B0  ,(ft)  with  0  <  B  <  1  and  *  r  ,  then 

P,d  p 

for  any  0  <  a  <  1  there  exists  a  piecewise  polynomial  tp  (x,y)  f  S^-’^,'''(ft”) , 
with  degree  vectors  P  and  £  in  which  p  -  =  q  .  =  1,  p,  ■  p.  .  =  p. , 

.  =  q.  .  =  p.  for  l<k<3,  15  IS  3,  l<j<  n+1,  and  p  =  max(2,[jy]) 

±fJ  ■L*J  J  J 

and  y  satisfying  (5.15)  such  that 


(5.11) 


|u-iji||  <  Ce“ 

Hx(ft) 


The  constants  C  and  b  are  independent  of  P_,  Q.  and  the  number  N  of 
degrees  of  freedom  of  the  space  S—  ^  (Q1!)  . 


Proof.  Applying  Lemma  5.1  there  is  (  S— ’-1*  (ftg.)  satisfying 


!|uJM|2  <c{  y  x  2(l-g)<  (pi,rSiJ)!  [hi 

u  ‘^*1  i  l  xj-l  \  (p.  .+S,  .-2) !  '  2 

H  (ft)  1<j— n+1  J  i,j  i,3 

l<i<3 


zzyrl^)  i,J  +3  2 

"•  H*1*!  ’Xj> 


(q.  .-t.  .)!  V.  .  2ti,j  ,,  ,,2  \  .  2(l-g)(n+l) 

+ _ 3-iJ.  1>2,.  Ht  +3  2  +a 

(q,  4+t4  .-2) !  I  2  J  „ti,3  ’  (a,  J 


i»J  3-0 


Letting  p  =  max(l,  — — )  we  have  using  (2.3) 


|u-4>!|2  —  c{  l  5 

H  (ft)  1  < j  5a+l 


LS 

j-1  1  j"2^  1  '  2  i.j 


+  — LilLLiJ-  (Pdj  *J((t  +l)*))2+a 

(q,  ,+t.  ,-2)!  I  2/  ((ti>:j+l).)) 


2  (1-6)  (n+1) 


Let  si,j  -  “j"!.! -  Vr  h.j  ’  “j’i.j  ■  Vr  0<  V*) 5 1  “ith  V“j 

being  determined  later.  Then 


>-12)  I!u-<PH\  5  C{  l  ( |  CF(pd,a. )  |"  ^  *P? 

HA(fl).  l<jSi+l  3  '  3  3 

+  |  F(pd,aj)|qJ^)  +a2a-B)(n+D} 


,1  ,1-a  /  ad..  r2a 

where  F(d1,ot)=-^2_  (  — J  • 

2  (1+a) 

and  F(pd,-j^j)  <  1. 


Function  F(.d^,a)  is  defined  on  IQ,00)  *  fO,ll 


Hence 


(5.13)  F 


min  *  «in  F(pdsa)  -  F<Pd,amin)  <  F(pd,  ^)  <  1,  amin 
a<E[0,l] 


Taking  p.  -  q.,  a  *  a.  “  max(-y-,  a  .  )  we.  have  by  (5.12) 
j  J  j  3  Pj  roin 


2/i+pV. 
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(5.14) 


HX(fl) 


<c[02(n+l)a-B)  +  £  x2(l-6)lF(pd>  _i_/dp6 

J  Pj  2 


+  I  x2(^-B)!f  p  pb} 

1+i  <1<n+1  J_1  111111  J 


l+j0<j^ri+l 


where  p  -f-M  +  1,  [a]  me 
J0  ^min 

max(2,[jvi])  for  j  ?  2  and 


means  the  integral  part  of  a.  We  choose  p.  = 


(5.15) 


>  2(1-B)lncr  _ 


Since  jQ  is 


bounded 


1  <  -  J-—  +  j 

J0  ~  u  a  .  m  ’ 
min 


(5.16) 


(tu-^l2  <  ca  2(n+1)  (1‘B) 

H  (fl) 


iPJ  6 


T  ,  p  uu.il  '  I  i 

1 


F(pd,-^)Pj  ,  ,p 


jQ+l<j<n+l  2  a2(1_6);i 


^  -  2(n+2)(l-0 

<  CO 


with  C  >0  independent  of  n.  Let  P  »  max(l,P).  The  number  of  degrees 
of  freedom  N  is 


N  s  3  +  3(2p+l)  +7  l  $ j+l)“  (n+2+=)  . 

j=3  y 

Hence 


IIM-II2! 

H  (SB 


-2  (1-g)  Ini- 
a 

<  Ce 


which  yields  (5.11)  with  b  = 


and  C  independent  of  N. 


Remark  1.  In  practice  p^  and  q..  are  often  selected  uniformly,  p^  * 
=  p  =  l  u(n+l)]  for  some  p  >  0.  Then  (5.11)  holds  and  (5.15)  is  not 
needed.  From  (5.16)  we  see  that  for  p^  =  p 


(5.17)  ||  2  <  Cc2(n+2)(1-B){l+|F  .  |Pp6  l  (i)2(1_S)j} 

TT^"  /a  \  lUill  1  _  m  _  ■  «  W 


1< jSn+1 


w<.«)(i*)  +  l 


2(1-6)11  3(S)  p3(f)3 

s  c!a  +  lF«J 


(f)2} 


<  Ce 


-bN 


-2  1 

/■i_  q)  3  13  1 

with  b  *  —  min(p  In  —  ,p  In  ■= - )  and  C  independent  of  N.  In 

3A  0  min 

(5.17)  the  first  term  is  the  error  on  the  element  containing  singularity, 


the  second  is  the  error  on  the  rest  of  the  elements. 


Let  us  discuss  now  the  optimal  choice  of  y .  First  we  drop  the 


term  p^  in  (5.17)  since  I ^minl  ^P^  -  CF1"  for  some  0  <  F  <  F^^  <  I, 

and  replace  and  y^^  by  N.  and  v  respectively.  Consider  the 

7  u  _2 


function  g(v)  =  g^(v)  +  §2^),  g^(\>)  =  ^^N0V  and  g2(v)  *  FV‘0^ 

with  0  <  F,  a<l-  It  can  be  readily  seen  that  g^v)  and  g2(v)  are 


increasing,  respectively  decreasing  functions  on  [  0 ,°°  )  »  and  g.(v)  ■  g0  (v) 

1/3 

for  v  *  Vq  *  2  ^ln^F^'  ~  with  Vq  independent  of  Nq.  g^(v)  and 


>2(v)  intersect  each  other  at  least  once  and  three  times  at  most,  (see  Fig. 


5.4).  Therefore  there  are  v^«  1  =  1,2,3  (resp.  i  =  1,2  and  i  *  1) 
such  that 

-2 

9  AM-ft'k 

(5.18) 


2(1-8)N0V 

2  4(1-0) _  In  Q 

Vi  *  N„v , 


F  °  iln  F 


since  g' (vQ)  >0  ,  g’(0)  <  0,  and  <  vQ. 


Figure  5.4 

Determination  of  the  Optimal  Value  of  V 


Now  we  can  speak  of  the  optimal  value  y  ■  (resp.  v.j)  in  the  sense 
of  minimizing  g(v)  at  (resp.  v.^  •  ^e  optimal  value  of  y  depend 

on  a.  F,  0,  and  N. 


For  the  asymptotic  analysis  we  write 


g(v) 


2(i-e)N0v‘ 


+  a 


In  F  „ 

i - viN. 

In  a  9 


If  0  <  v  <  Vq,  g(v)  z  a"1"  a  for  large  N^,  which  is  decreasing  with  v. 

2(1-0)N  v"2 

If  v  >  Vq,  g(v)  ~  a  ,  which  increases  with  v.  Hence  v0  is 

the  asymptotic  value  of  the  optimal  V.  For  large  N  equation  (5.18) 
has  only  one  solution,  and  will  approach  to  from  left  and 

right  respectively. 


Remark  2. 


a  .  and  F  .  are  the  function  of  d,  *  pd,  and  F  , 
mxn  min  1  min 


F(dl-“nin  <V>- 

We  see  that 


Fmin(dl>  '  (Fd1  +  «in(V  V 


>0, 


and  F  ^  is  increasing  function  of  d^.  If  p  and  d  small,  we  can 
choose  small  y ;  hence  we  can  reduce  the  number  of  degrees  of  freedom  of 
sP.2,.1 

(ft") .  The  value  d  is  related  to  the  function  u(x,y)  and  P  is 
related  to  the  mesh. 


Remark  3.  The  mesh  ft*"  with  the  irregular  points  can  be  used  analogouslv. 

-  a 

The  results  are  then  quite  the  same  as  for  the  triangular  mesh  we.  analvzed 
above. 


Remark  4.  The  technique  we  used  for  the  adjustment  of  two  adjacent  elements 
can  be  directly  applied  to  satisfy  homogeneous  (or  polynomial)  boundary  condition 


So  far  the  singularity  of  the  solution  occured  in  one  vertex  only. 

The  general  case  can  be  easily  treated. 

Let  u  F  Bg(fi)  with  8  =  (8^62*63.64),  0  <  ^  <  1,  1  <  i  <  4  and 
ft  *  (0,2)  x  (0,2).  Divide  fl  into  four  subdomains  mft,  1  5  m  5  4  such 
that  on  each  mft  the  function  u  has  only  one  singular  point,  and  apply 
Theorem  5.1  on  1  <  m  <  4.  Let  be  the  union  of  the  geometric 

meshes  Mft£  on  ,  1  <  m  <  4.  Let  P  =  (1P,2P,3P,4P)  and  =■  (1£, 2£, 30 , 4£) 

be  the  union  of  mP  and  m£,  1  5  m  <  4,  each  ^P  and  be  defined  as 

before.  Thus  for  the  mesh  and  degree  vectors  P_  and  (£  (see  Fig.  5.5) 

we  have 

Theorem  5.2.  Let  ft  =  (0,2)  x  (0,2),  ft^,  IP,  (£  be  defined  as  above. 

If  uFB2  d(ft)  with  6  =  (8^62,83,64),  0  <  3±  <  1,  1  <  i  <  4,  then  for 

any  0  <  a  <  \ ,  there  exists  a  piecewise  polynomial  \p  (x,y)  (  S^’-^’*(ft”) 

such  that 

2  h„l/3 

l|u-^||  x  <  Ce  bN 

ir(fl) 

—  1/3  U/3 

with  constant  C  independent  of  N  and  b  *  2(l-6)ln  —  1  I  where 

0  »28pz  ' 

8  *  max  g . ,  U  -  max(l,y)  and  ]i  satisfies  (5.15).  0 

lSi<4 

We  shall  now  briefly  address  the  geometric  mesh  on  a  parallelogram 
domain  ft,  which  is  the  image  of  geometric  mesh  on  square  under  a  linear 
mapping.  Let  M  be  a  mapping  of  standard  square  S  =  (0,1)  x  (0,1)  onto  ft: 


8  =  (31,g2,B3,g4) 


(5.18) 


x  =  a£  +  bn  cos  co 


^  y  ■  bn  sin  to 

with  Jacobian  =  a  b  sin  u).  The  inverse  of  H  is  given  by 


(5.19) 


£  =  (x,y)  =  (x  -  y  cot  w)/a 


n  “  \l/0  (x,y)  =  y 


sm 


aj/b 


where  u)  is  the  interior  angle  of  ft  and  a,  b  are  the  lengths  of  the 
sides  of  ft  shown  in  Fig.  5.6.  Let  ft^  *  M(Sn) ,  where  S^.  is  the 
geometric  mesh  of  square  S  (see  Fig.  5.2). 

Let  S^og)  be  the  subspace  of  continuous  and  piecewise  polynomials 

4)(x,y)  (  (ft)  such  that 

d>|  ,  7  C.  A®  for  1  <  k  <  3,  2  5  i  <  3,  2  <  j  <  n  +  1, 

k  L  £,mrlr2 

“J  *  ^ 

» J  Q<£<p .  . 0«m<q .  . 


and 


<l>l 


ft 


i,j  0<£Spi^,  0<m<q_ 


for 


1,  1  5  j  5  n  +  1  and 


i  =  2,  3, 
j  =  2. 


The  number  of  degrees  of  freedom  of 


s£*a,1(s“) 


is  defined  as  the  number 
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of  degrees  of  freedom  of  S— and  denoted  by  N.  We  have  the 
following  theorem 


Theorem  5.3.  Let  D  be  a  parallelogram  shown  in  Fig.  5.6  with  interior 
angle  to,  0  <  to  <  -g-.  Then  for  any  u  $  B2  ,  (fl)  with  =  r^,  0  <  g  <  1, 

^  P,h  p 

and  0  <  a<l,  there  exists  t//(x,y)  f  with  the  degree  vector  P 

u 

and  £  in  which  p  =  q  =  1,  p  «  q  •  p^  =  -  maac(2,  [jul), 

x,  x  x,x  x»j  x,j  i,j  x*j 

j  -  2  and  p  satisfying  (5.15),  with  3  replaced  by  3  =  d/2  max(a,b)/ 

min(a,b)  such  that 


lu-.ll2,  5  Ce-bN 
H  (ft) 


with  b  »  2(1-6) 


&T  *  i 


,  y  =  max(l,y)  and  constant  C  independent  of  N.  □ 


b  sin  cj 


^3,4 


a  b  sin  w 


a2  b  sin  w 
o-3  b  sin  w 


AY”  / 
MaY®. 

r*h, 2/&U2 


0  a<r3  oar2  a  c 


Figure  5.6 


Geometric  Mesh  on  a  Parallelogram 


6.  GENERALIZED  GEOMETRIC  MESH  AND  ITS  APPROXIMATION  PROPERTIES 


In  this  chapter  the  generalized  geometric  mesh  with  triangular  and  quadri 

lateral,  curvilinear  triangular  and  curvilinear  quadrilateral  elements  will 

be  employed  on  polygonal  domain  £2  contained  in  the  unit  disk  centered  at 

the  origin  which  coincides  with  one  of  the  vertices  of  £2.  We  need  to 

redefine  mesh  £2n  and  space  s— ’^(£2n). 

c  a 

Let  £2°  =  {£2,  j  =  l,2,...,n+l,  i  =  1,2, . . .  ,I(j) }  be  a  partition  of 
*  J 

£2  satisfying  the  following  conditions: 


1)  0i  's  are  quadrilateral  or  triangles  (curvilinear  quadrilaterals 
or  triangles).  The  intersection  of  any  two  £2  's  is  one  common  vertex,  or 
one  entire  common  side  or  is  empty; 

2)  Let  h.  .  and  h.  .  be  the  maximum  length  and  the  minimum  length 

1 » J  1 » 1 

of  sides  of  £2  .  and  for  all  i,j  there  is  a  constant  X  such  that 

1 » J 


(6.1) 


h.  Jh,  .  <  x; 

“i»J 


3)  Let  M  *  (M.  ,  IS  i5  I(j),  1  S  j  <  n  +  1)  in  which  M,  .  is  a 

1 » J  1  *  J 

one-to-one  mapping  of  standard  square  S  =  (0,1)  x  (0,1)  (resp.  standard 


triangle  T={(£,n)|0<n<l-£,  0  <  £  <  1}  onto  £2.  ..  If  P.  and 

1  ^  ^ 

y  denote  the  vertices  and  sides  of  £2.  . ,  then  M  (P  )  and  M-  ,  (Y.) 

X*  i,  J  i.j  x,  ^ » J 


are  the  vertices  and  sides  of  S,  1  S  l  S  4.  Moreover,  if  M.  .  and 

i.J 


M  ,  map  standard  square  onto  two  elements  £2.  .  and  £2.  ,  with  common 
m,k  r  n  i,j  2,k 


side  y  =  AiAo  then  for  any  A  ?  y,  dist(M.  .(A),  M.  . (A«))  =  dist(M  .  (A)  , 
■L  £  i  >  J  ~  ni,K 


-1 


Mm  k^o^’  1-^-2.  It  is  assumed  that  the  mapping  can  be  written  in  the 


form 


x  =  Xi  j(5,n) 


(6.2) 


on  S(resp.  T) 


y  -  Y4  ,(5.n) 


with  X.  .  and  Y.  .  being  smooth  functions  on  S, 
i » 3  i  »3 

(6.3)  |Dax|,|Day|  <  Ch.  |a|  <  2 

1 » 3 

and 

(6.4)  C.h2  .  <  J.  .  <  C.h2  . 

1—l,2  1,3  2  i,3 

where  J_^  is  the  Jacobian  and  the  constants  C,  and  are  independent 

of  i,j . 

The  mesh  Qn  (0<cf<1)  is  called  the  geometric  mesh  if  in  addition  the 
a 

following  condition  is  satisfied: 

4)  If  d.  .  denote  the  distance  between  the  origin  and  Q.  ,  then 

(6.5)  an+2~j  <  d  <a  n+1~\  1  <  j  S  n+1,  1  <  i  <  I(j), 

x  y  J 


(6.6)  dA  1  -  0,  1  <  i  <  I(j), 


(6.7)  k  d.  .  2  h.  .  5  h.  .5  k  d.  .  for  1  <  i  <  I(j),  1  <  j  5  n+1, 

1  2  i,j 

(6.8)  ^  ~  — i  i  -  i  -  for  1  5  i  5  1(1) 

where  K^,  1  5  i  <  4  are  the  positive  constants  independent  of  i  and  j. 

Let  P  =  (p.  1  Si  2 1  ( j )  ,  l<j<n+l)  and  Q  =  (q  .,  l<i<I(j), 

1  »3  I » J 

lSj  Sn+1)  be  the  degree  vectors  with  integers  p  ,  q  >0. 

p  i’j  1,2 

We  define  the  subspace  S~ ’^(Qn)  as  in  [9] 


S-’%f)  =  { cf)  |  <J)(x,y) 


<f>.  .  (M.  .  (x ,y) )  for  x  €  Q.  ., 
i » J  ,  J 


<j>.  ,  (£,n) 

j-.j  j-jj  j->j  i>j 

is  the  polynomial  of  degree  Sp.  .  in  ?  and  of  degree 

i  *1 

5q.  .  in  n  on  S} 

1 » 1 

and  S-’^,;L(fin)  =  S-’-^(Qn)  n  H1^),  (£)n)  =  S -’Q*1  (tia)  fl  H1  (u) 

O’  O’  o  o 

where  H^(£2)  =  (u  |  u  €  H^ft),  u|  ^  -  0}  . 

Remark  1 .  Denoting  U  =  u(M.  .  (C,0)),  (6.3)  implies  that  for  |a|  =  1 

1 » J 


(6.9) 


|Dau|  £  ChT1  (|U  |  +  [U  ()  on  S, 

1  >  j  s  n 


and  for  a  =2 


(6.10) 


|d“u|  =  Ch-^J  <|uKl  +  |unnl  +  |«5nl  +  |Un|  +  |u,|). 


Remark  2.  If  £2  is  a  triangle  ,(6.3)  and  (6.4)  are  equivalent  to  the  well 
-  i,j 

known  angle  condition: 


(6.11)  0  <  Ug  5  u  5  ir 


“0  "  * 


for  all  interior  angles  io's 


of  £2.  .  . 


Remark  3.  The  geometric  mesh  £2n 

-  a 

2  g 

functions  u  £  B“  (£2) ,  $  =  r  . 

D  ,d  p 

for  the  approximation  of  functions 


was  designed  for  the  approximation  of  the 
In  an  obvious  way  the  mesh  can  be  designed 


u  €  (ft) , 

p  ,d 


M 


Lemma  6.1.  Let  £2  be  a  curvilinear  quadrilateral  (resp.  triangle),  h  be 
the  length  of  the  longest  arc,  and  let  M  be  a  one-to-one  mapping  of  standard 
square  S  =  (0,1)  x  (0,1)  (resp.  standard  triangle  T  ®{  (  5,0  )  !  0  <£  <n  ,  0  <£<!}) 
onto  £2  given  by 
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(6.12) 


x  =  <f>(C,n) 


y  =  ifKS.n)* 


on  S  (resp.  T) 


Assume  that  for  any  a,  jal  Si 


(6.13) 


|Da(p|,  |Da^|  <  Alalh. 


Then  for  U(£,n)  =  u(cj)(^,n) »  iKC.n))  and  (a |  =  k 


(6.14)  ||OaU||  Q  0  <  (3A)  £  max  { [[D^u) 

H  (S) (resp.  H  (T))  2=1  |y|=2 


,  h£(k-l)! 

0  *  (o-p ! 

H  (Q)  ' 


where  y  is  the  multi-index  (y^,^) 


Proof .  Let  us  first  show  that  for  | ct  |  =  k  >  1 


(6.15)  D0!!  -  l  l  (D^u)p^k) 


*  1  y^s* 


where  Sx  is  a  set  of  pairs  y  =  (y^,y^)  y  ,y  >0  integral,  y  +  y^ 


such  that 


s  =  (1,0)  u  (0,1) 


k+1  k  k 

(6.16)  S  =  { (y  ,y  )  j  (y  -l,y  )  (  S  ,  (y  ,y  -1)  €  S  } 

PIP  q  PI 


S  consists  of  pairs  which  are  in  general  repeating.  For  example 


=  {(2,0),  (1,1),  (1,1),  (0,2)} 

S3  =  {(3,0),  (2,1),  (2,1),  (1,2),  (2,1),  (1,2),  (1,2),  (0,3)}. 


The  values  (y  ,y  )  are  the  numerical  values  of  the  a  pair  y  and  the 

P  q 

different  pairs  can  have  the  same  numerical  value.  Further  we  will 


(not  only  its  numerical 


denote  j  y [  *  +  u  . 


(k)  — 

The  functions  p—  depend  on  the  pair  u 
^T,  a 

value).  For  example  for  Z  =  2,  a  =(1,1)  we  have  in  (6.15) 

y  (dHi)o^  =  u  d>d)_  +  u  A  ip  +  u  +  u  Wj 

.  o  u  n  xx  S  S  xyvC  n  xyvcvn  yywcvn 


uss2  p’a 

(2) 

and  hence  the  function  p  associated  to  the  two  different  pairs  with 

U»a 

the  same  numerical  values  (1,1)  are  different.  In  (6.15)  the  general  form  of 

(k)  . 
p_  is 

U»a 

^  V .  S2  K . 

(6.17)  p_^k)  =  l  _  a(]j,v.  ,ic.  ,a)  TT  D^4>7TD^ 

U.ct  s^s^lhl  J  2  j=l  j=l 

lVjl+lKjl=k 


with  a(y,v  ,K.,a)  =0  or  1.  (We  used  in  (6.17)  the  notation  "]"T=  1) . 
J  J  i 

(6.15)  can  be  readily  proven  by  induction.  Let  $(k,£).k  >  Z  be  the 

number  of  (additional)  terms  in  (6.17).  It  is  immediate  that 


(6.18)  4(k,k)  =  1 
and 

(6.19)  $(k+l,£)  =  Z<P(k,Z)  +  2$(k,£-l) 

and  (6.18)  (6.19)  gives 


4(kt4)  <  3k£k"S'. 


Hence  coming  back  to  (6.17)  we  get  by  (6.13) 


and  hence 


|Dauil  0  <  (3A)k  I  max  {1|d“u||  q  >hV~S 

a  (S)  A-l  |u|-z  H  (ft) 


-  (3A)k  l  max  {||Dau||  Q  >h£ 


1=1  |y|-A 


Hu(ft) 


2  ft 

Theorem  6.1.  Let  u  €  B„  j(ft)  with  $„  =  r  ,  0  <  Q  <  1,  and  Q  be  a 

p,d  p 


polygonal  domain  contained  in  a  unit  disc  and  with  a  vertex  at  the  origin 


Let  ft”  =  {ft.  .  ,  I2i51(j),  1  <  j  <n+l}  be  a  curvilinear  quadrilateral 
o  1 ,  j 


geometric  mesh  satisfying  conditions  l)-4).  Assume  that  mapping  M. 

1  >3 


S  =  (0,1)  x  (0,1)  onto  ft.  . 

J 


f  x  38  <f>.  .(5,0) 

1  9  J 


(6.20) 


in  S 


y  -  ^^(5,0) 


is  such  that 


(6.21) 


1  Da4) .  ,i,|DV  A-  Akh.  *  for  any  |a|  -  k>  1,  1  <  i  <  I(j)  ,  1  <  j  <  n+1 
1  *  J  -LjJ  1  *  J 


where  A  is  independentof  i,  j,  and  k.  Then  there  exists  a  function 

cj>(x,y)  €  S— ’^-’^(ft11)  with  degree  vectors  and  Q  in  which  p.  .  =  q 

o’  i ,  j  i  i 


Pj  =  max(2,  [ j M ] )  ,  j  >  2,  and  with  y  satisfying  (5.15),  Pi  i  *  ^ 


such  that 

(6.22) 


|u-4>||  i  5  Ce 

H  (ft) 


-bN 


1/3 


with  C  and  b  independent  of  N. 


Proof.  For  each  (i,j),  U.  .(5,n)  ■  u(4>.  .(£,n))  is  defined 

-  i»J  i,j  i,j 

on  5  =  (0,1)  x  (0,1)  and  by  (6.21)  and  Lemma  6.1  the  inequality  (6.14)  holds 

Applying  Lemma  4.3  there  is  a  polynomial  $ .  .  (£,n)  of  degree  p.  in  £  and 

3 

H  such  that  for  the  integer  t.,  1  <  t.  <  p.  and  0  5  m  <  2.  i  >  2 

3  3  3  ~ 


m  (p,-t,)! 

ID  (U.  )||  n  <  C  —  xr~Zo~9~T 

'  X»J  H°(S)  (Pj+tj+2-2m) 


t.+l+k 


. ,  t.+l+k 

2(t.+x)  9  a  3  T 


r<|)  l  (II  t4.>¥0  (i)' 


k'°  H3  H  (S) 


by  (6.14) 


3  J  U.  .  0  .  2k 

+  'I — r+iJ||2o  <|>  > 

3*n  3  H  (s) 


2(t.+l)  t.+3 


.  c  (W!  -3A.  J  T  „2*  2„  .  .2 

-  (pj+tj+2-2m)!(  2>  ^  i.jW  I|D  “”H0 

i » J 


by  (6.7)  and  because  d  <  1 

i.J 


(Pi-tj)!ti  3A  2(tJ+1) 

— _J _ 3 _ ] _ 

(p.+t.+2-2m) r  2 ’ 

3  3 


d2(2-e)Cf3 

i,J  W  *2  Wh2,2([!)- 


Letting  (j).  .(x,y)  =  .  (M  (x,y))  for  j  >  2,  we  have  by  the  scaling 

it*!  1»J  1 


argument 


|Dm(u-4>:  ./II2 


10  H  (fl.  .) 

i.J 


C(p  -t .)  !t , 


(Pj+tj+2-2m) 


-AZti  d2(2-m-«  f  ’  )2 ||U||; 

i  2  di,j  /  2  Mi-1)!  }  ,lU|1  £,2 


Hg  («) 


Let  H0’  (ft)  (k-:2)  be  the  weighted  Sobolev  space  with  norm 


<sa 


■ll-Mf  9 

C’(ft) 

P 


v  2&,(k-l)f.2  ,,  2 

ih  *  «'»!  lluV-2(n> 

P 


Since  u  €  ,(fl)  then  for  k  >  2 

P>“ 


MLk  2  i  Cd1'k"2)(k-l)! 

He‘  <n) 

A.  AlfiO  O 

with  d  =  max(l,K  d).  Let  T  be  the  operator  H  ’  (Q) 

p 

1  <  k  <  p,,  0  <  m  <  2 ,  Tu  =*  u  -  4>  ■  ,  then 

J  i,j 


H  (ft.  .)  for 


Hk+3’2(fi),Hm(fl  .)  ^(pj 
i » J 


<  /  (pi~k)!  /3A\2k  2(2-m-3)\1/2 

~  L\(p,+k+2-2m)!\  2)  di,j  )  * 


t  +3  2 

Let  H  3  (ft)  *  (Hk+2,2(fi)  ,Hk+3,2(Q))  be  the  interpolation  space 

P  P  P  o ,°° 

by  the  K-method  [10]  for  t.  =*  k  -  1  +  0.  Then  T  is  linear  operator: 

V3’2  m  3 

J  (ft)  -*■  H  (fl.  .)  with  norm 
P  p 


t  +3,2 

i  J  um 


(Q)  ,H  ) 


r(pi~t1+1)  f3A\2tj  2(2-m-0) 
r(p.+t.+3-2m)  \2  >  ai,j 


(see  the  proof  of  Theorem  4.1).  Thus  for  a  real  number  t^ ,  1  S  t  5  p 

(6.23)  0D“(u-*  „20  *  . 


By  (2.3) 


(6.24) 


U  *-  i  o  o 


1/2  Ci+^ 

<  ct;^  d  2  r(t  +2) . 
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For  j  *  1,  <J>  is  always  chosen  as  bilinear  interpolation  of  U,  ,  at 
i » i  i  9 

vertices  of  $2.  .  .  Letting  <p.  ■,  <J> .  ,  0* .  1(x,y))  we  have  by  Lemma  3.6 

i,I  I,*  i,-*- 

and  the  standard  scaling  argument 


K  J  i  *  cj|u.  -*  ||  1 

1,1  H  (A.  .)  1,1  1  H  (S) 

^  *  J- 

-  c|Bi  r¥ii  2  2 

1,1  »8  <S) 


where  $  = 


'vg  'V  — 1 

r  ,  r  is  the  distance  between  point  (C»h)  €  S  and  M  (0,0), 

1,1 


(6.25) 


f|u~0±  ,11  1 

'  *  «1,1> 


S  2  2 

B  i,l 


To  achieve  continuity  of  the  polynomials  on  two  adjacent  elements  we  use  the 
same  procedure  as  we  used  in  the  previous  chapter.  We  construct  in  this  way 
the  function  <j>  €  S— ’^^(n^)  such  that 


(6.26) 


Mil*!  hi(i"e)  l«l^2  2 

H  «!>  i»l  1,1  tC 


I,. .(f) 


2t.  r(p  -t  +1)  2(l-R^  2 

J  d"(1  e)iiuir 


^(3)' 2 '  rwi)  1-j  V/3’2, 

l<jSn+l  Hg  1 


S  C02(1_0)(n+2){I(1)  +  l  I(j) 

j-2 


[3Ad,q) |  j  6. 

2(1-6) j  V 


where  t.  *  a.p.,  and  F(d, ,a) 
J  j  J  1 


(l+a)1+a  '  2  ' 


It  has  been  proved  that 


min  F(d,  ,ot)  -  F(d., —  -  -) 


a 


V 


F  <  1. 
min 


4+d! 


It  can  also  be  shown  by  (6.5),  (6.7)  and  (6.8)  that  X ( j )  J  K,  for 
1  <  j  <  n+1  and  some  K  >  0  independent  of  N.  Letting 


a.  =  max  ^  — 

3  Pj  ( 3  Ad ) 2 


•},  Pj  =  ®ax(2,  [jy]) .  j  >  2  with  y  satisfying 


(5.15)  we  have  from  (6.26) 


|u-<{>j|‘\  <  CKa 

H  (fl) 


2(1-6) (n+2) 


<  Ce 


-2(l-6)(-T2)1/3  ln(~)  *Ni/3 
K U  c 


with  C  independent  of  N  and  y  =  max(l,y). 


2  B 

Theorem  6.2.  Let  u  €  ,  ((2)  with  =  r  ,  0  <  6  <  1,  and  S2  be  a 

p,d  6 


polygonal  domain  (2  contained  in  a  unit  disc  and  with  a  vertex  at  origin. 


Let  (2^  =  {(2^  i  1  S  i  S  I(j)i  1  £  j  5  n+1}  by  a  curvilinear  triangular 


geoemtric  mesh  satisfying  conditions  l)-4)  in  Section  4.4.  Assume  that  the 
mapping  M.  .  of  T  =  {(£,q)  |0<n»  0<£,  £+n<l}  onto  (2. 


x  =  <|)i^(C,n) 


in  T 


y  =  4>,  , 

»  j 


is  such  that 


|Da<}>.  .  |  ,  |D°fc  .  |  £  A^h.  ..  for  any  |a|  =  k  >  1,  1  <i  <I(j),  1  < j  2n  +1, 


where  A  is  independent  of  i,  j  and  k.  In  addition  assume  that  . 


can  be  extended  to  standard  square  S  »(0,1)  x  (0,1),  M  (S)  c  n  preserving 

^  >  J 


the  properties  of  M 


on  T  and  that 


for  all  i  and  j 


d  .  *  dist(M  . (S) ,0)  >  <  d  . 

i,3  0  i ,1 


with  constant  >  0.  Then  there  exists  <p  (x,y)  €  (S~ ®  S— 

with  degree  vectors  P  and  Q  in  which  p.  .  =  q.  .  ■  p.  =  p.  *  max(2,[jy]) 

1  »3  i* J  1  J 

for  j  >  2  and  y  satisfying  (5.15),  p  =  q.  =1,  and  with  degree 

i ,  i  x ,  i 

vectors  ?  and  ^  in  which  ’p  .  =  1,  q.  .  =  2p.-l  such  that 

i»3  J 

1/3 


l«-<HI  ! 


<  Ce 


-bN 


H  (fl) 


with  C  and  b  independent  of  N. 


Proof .  The  proof  is  analogous  to  the  proof  of  Theorem  6.1. 

Remark  1.  We  assumed  in  Theorem  6.2  that  the  mapping  can  be  extended.  In 

contrast  to  the  mapping  of  T  where  one-to-one  mapping  is  assumed  together 

with  condition  4)  we  do  not  need  these  assumptions  to  be  satisfied  on  the 

entire  S,  we  need  only  that  M.  , (S)  c  ft. 

1  >3 

Remark  2.  Assume  that  ^  is  a  simple  parallelogram  or  triangle  only  and 
2 

assume  that  u  ^  Bq  d  ^Ut  suc^  t^iat  domain  of  analyticity  can  be 
extended  into  a  small  neighborhood  of  ft.  If  ft  is  a  parallelogram 
and  the  mesh  consists  of  one  element  only  then  the  rate  of  the  p-version  is 
exponential.  It  follows  immediately  from  our  analysis.  If  ft  is  a  triangle 
and  the  domain  consists  of  one  (triangle)  element  only  then  we  cannot 
conclude  directly  from  our  theory  that  the  rate  of  convergence  is  exponential. 

Remark  3.  In  Theorems 6.1  and  6.2  we  assumed  that  the  domain  is  a  polygon. 
Obviously  this  assumption  is  not  essential. 


Remark  4.  As  in  the  previous  chapters  we  can  consider  in  Theorems  6.1  and 
6.2  uniform  degree  of  elements  and  get  the  exponential  rate  of  convergence  too. 


Remark  5.  The  quadrilateral  is  the  special  case  of  a  curvilinear  quadrilateral. 
Suppose  that  (x_^,y^)  are  the  vertices  of  Q.  The  mapping  of  S  =  (0,1)  x  (0,1) 
onto  is 

j  x  =  iK£,n)  =  xi  +  (x2-xi^  +  (x4-x1)ti  +  (x1-x2+x3-x^)5n 
y  =  *KS,n)  =  y1  +  (y2"yi^  +  (y4~yi)n  +  ^i~y2+y3~y^^Tl- 
Obviously  .  2  1)  holds. 

Remark  6 .  In  the  mesh  shown  in  Fig.  6.1  the  curvilinear  triangular  and 
curvilinear  quadrilateral  elements  are  combined.  This  kind  of  mesh  is 
important  in  practice. 


Figure  6 . 1 

Curved  Geometric  Mesh 


Figure  6.2 

Mapping  of  Standard  Triangle 


The  mapping  of  S  =  (0,1)  x  ( 
element  is 

x  =  4>(5,n) 

^  y  =  iKS.n) 

with  h  =  qir/4  and  (6.21)  is  satisfied. 

The  mapping  of  T  onto  curvilinear  triangle  (element  2)  can  be  defined 
in  various  ways .  Let  us  mention  one 

^  x  =  <(>(£>n)  =  q2(n  cos  (cos  -  q  cos  ^-)) 

< 

=  'K€,n)  *  q2(n  sin  ^  T  -£§^  (sin  ~  n  sin  ^)). 

Observing  that 

/  nn  i 

(cos  ^  -  n  cos  J)  “ 

and 

,  .  TTTl  .  TT.  1 

(sm  —  "  h  sm  — 

2 

is  analytic  at  q  =  1  and  we  easily  see  that  with  q  *  h  (6.21)  is 
satisfied.  If  is  also  easy  to  see  that  M(S)  c  ft  (see  Fig.  6.2). 

Remark  7.  Consider  an  element  ft  with  vertices  (x_^,y^),  1  -  i  -  4  and 

curvilinear  edges  y^  (see  Fig.  6.3),  y^'s  are  described  by 


,1)  onto  the  curvilinear  quadrilateral 

*  7  h(5(l-q)+q)  cos  ~r 

<1  4 

=  ^  h(£(l-q+q)  sin 

TT  4 


1  ~ 

xi ' 

g.(?) 

/ 

\ 

for 

0  s  5  <  1,  i  -  1,3 

U  ■ 

fi(0 

p  ■ 

gi(o) 

for 

o<n<i,  i  =  2,4 

Cyi  " 

fi(n) 

gi  and  fi. 

1  S  i  S  4 

are  infinitely  differentiable,  and 

1  (k>l 
|8±  1  * 

|£<k)|  <  h 

Ir 

A  ,  for  IS  iS  4,  any  k  >  0 

and 

IVgi+2>’  lfi"fi+2l  -  Ah  for  1  =  1»2 

with  A  >  1  independent  of  k,  h  and  h  being  the  arc  length  of  the 
longest  curvilinear  edges  of  Then  the  mapping  of  S  =  (0,1)  x  (0,1) 

onto  u  by  a  blending  function  is  constructed  as  follows: 

!  x  =  <?>(£, n)  *  x^Od-n)  +  x2(nK  +  x3(On  +  *4(n)  (i-£) 

-x1(i-C) (l-n)  -  x2£(i-n)  -  x3£n  -  x4n(i-0 

/ 

]  y  -  tKC.n)  -  y^Dd-n)  +  y2<n)  +  y3(S)n  +  y4(n)  d-5) 


1 


-y.  d-O  (l-n)  -  y-Cd-n)  -  -  y,n(i-D 
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3 


S 


11 

0 


21 

I 


Figure  6.3 

Mapping  of  Standard  Square 

Then  it  is  easy  to  see  that 

$  =  U-n)x{(£)  +  nx^(c)  +x2(n>  -  *4(n)  +  (x1-x2)d-n)  +  (vx3)n’ 

|4>5I  -  5  Ah; 

<|>  z  -  (i-n)x^}(5)  +  nx^CC).  -  -xJA)(0  +  x^d) 

|4>  »l  .  |<fr  JJ,  i  -  2AJ°h ,  for  l  >  2, 

5 

^  =  -x^(0  +  x3(0  +  (X2“X1)  +  (x4_Xl)’  -  4Ah; 

eh  „  =  0  for  m  >  1,  l  >  1  or  l  >  2,  m  >  1. 

C  n 


Similarly  we  can  see  that 


<j>  I  ,  <j>  2  2A  h  for  m  >  2, 

m  -  m 

n 

U  I  <  5 Ah. 

I  I 

These  inequalities  are  also  true  for  'p.  Hence  (6.21)  holds. 

Remark  8.  In  the  theorems  above  we  assumed  that  the  singularities  are 

i 

A.  are  the  vertices  of  (2.  Assume  now  that  A.  are  located  outside  of  £2. 
i  i 

*“cL  I/O 

Then  it  is  possible  to  show  that  ||e||  5  Ce“  N  '  .  This  rate  of  convergence  is 
achieved  by  the  p-version  when  the  size  of  the  minimal  element  of  the  optimal 
mesh  is  not  going  to  zero  as  p  -*■  00 . 


located  in  the  vertices  of  the  domain 


i.e.,  $  (x)  =  7T!x  -A 

3  Z-l 


Bi 

when 


7 .  NUMERICAL  RESULTS 


In  this  chapcer  we  will  discuss  the  numerical  results  for  the  solution 
of  a  model  plane  elasticity  problem.  We  selected  the  model  of  a  cracked  panel 
loaded  by  such  tractions  that  the  exact  solution  is  the  first  and  second  mode 
of  the  stress  intensity  factor  solution.  This  problem  was  selected  because 
it  characterizes  the  difficulties  of  the  usual  engineering  problems.  We  will 
compare  the  performance  of  the  h,  p  and  h-p  versions  of  the  finite 
element  method  by  focussing  on  the  accuracy  measured  in  the  energy  norm. 

The  purpose  of  the  numerical  computation  is  the  following: 

1)  Our  estimates  are  upper  estimates  which  have  asymptotic  character. 

It  is  important  to  see  the  numerical  behavior  of  the  error,  its  asymptotic 
range,  the  size  of  the  constants  characterizing  the  error  and  the  maximal 
accuracy  which  is  practically  achievable. 

2)  The  h-p  version  is  characterized  by  the  mesh  factor  a  of  the 
geometric  mesh  and  the  degree  factor  y  governing  mesh  size  and  the  growth 
of  the  element  degrees,  respectively.  Numerical  results  will  show  the 
sensitivity  of  the  accuracy  on  o  and  y,  the  values  of  the  optimal  factors 
a  and  y  leading  to  the  highest  accuracy. 

3)  Our  theory  does  not  allow  us  to  distinguish  between  performance  of 
elements  of  various  type,  the  curvilinear  and  rectangular  (resp.  triangular). 
The  computer  time  is  smaller  for  rectangular  or  triangular  elements  than  for 
curvilinear  elements.  Therefore  the  question  arises  whether  curvilinear 
elements  in  the  neighborhood  of  tiie  c rack  can  improve  the  accuracy  because  the 
singularity  has  radial  character. 

4)  It  is  known  that  the  p-version  is  very  insensitive  to  the  size  of  the 
Poisson  ratio  (nearly  incompressibility)  (see  [8  ],  [22]).  Numerical 
solution  will  show  the  effect  of  the  Poisson  ratio  on  the  accuracy  of  the 
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solution  in  the  case  of  the  h-p  version, 

5)  The  standard  finite  element  programs  are  based  on  the  h-version 
with  low  order  elements.  Several  hundred  finite  element  programs  have 
already  been  developed.  One  of  the  most  popular  packages  is  NASTRAN 
in  its  various  versions,  for  example  MSC/NASTRAN.  Other  widely  used  packages 
are  for  example,  ADINA,  ANSYS,  STRUDL,  GIFTS,  PAFEC ,  etc. 

There  is  only  one  code  based  on  the  p  and  the  h-p  versions,  the 
code  PROBE.  The  architecture  of  this  program  is  different  from  the  above 
mentioned.  It  utilizes  the  hierarchic  type  of  elements,  computes  simul¬ 
taneously  solutions  of  different  degree  elements,  etc.  We  will  make  some 
(crude)  comparison  of  computation  by  PROBE  with  the  fr-version. 

The  theory  and  computation  addresses  only  the  performance  with  respect 
to  the  error  measured  in  the  energy  norm.  Although  in  practice  other 
measures  are  essential,  the  energy  norm  performance  is  obviously  the  starting 
point  of  main  theoretical  importance. 

The  computation  of  the  h-p  version  has  been  done  by  the  program  PROBE  of 
NOETIC  TECHNOLOGIES,  Inc.,  St.  Louis.  (See  [21]). 

The  h-version  computation  has  been  done  by  the  adaptive  program  FEARS 
developed  at  the  University  of  Maryland.  (See  [14]  [19]). 

We  shall  consider  the  plane  strain  problem  of  two  dimensional  elasticity 
(homogeneous, isotropic  material)  with  E  and  v  denoting  the  Young's 
modulus  of  elasticity  and  Poisson  ratio  respectively  (E>0,  0-v<.5).  The 
domain  D  under  consideration  is  a  square  panel  with  a  crack  as  shown  in 
Fig.  7.1. 


Figure  7.1 
Cracked  Panel 


Let  U  »  (u,v)  be  the  displacement  vector,  t  = 


/Tll  X12\ 

~  \T21  T22 J 


be  the 


stress  tensor  and  T  be  the  traction  vector.  The  displacement  vector  U 
satisfies  the  Lame-Navier  equations 


(7.1) 


■<x+“>  vA  *  g> 


....  3  ,3u  3v. 

■(X+|J)  + 


in  D 


and  the  boundary  condition 


(7.2) 


T'3D  *  T  *  ‘ 3D 


where  n  is  the  unit  outside  normal  to  the  boundary  3D.  X  = 


(1+v) (l-2v) 


and  p  *  2 (1+v)  are  t*ie  ^ame  coefficients.  The  hilinear  form  associated 


with  (7.1)  and  (7.2)  is 
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(7.3)  B(U,W) 


an  ^wi  3W9  \  a..  a,.  ^w’  3w, 

r,.  /■£H _ _  .  3v  _ 2\  t  A/£i±  ,  °vn  , _ i.  ,  _ £n 

3x  3x  3y  3y  2  3x  3y^  ^3x  3y  J 


,  n  /3v  3uw^W2  ^Wlu,  , 

+  2  (57+  3?)(sr+  3F>>d,tdy 


with  U  »  (u,v),  M  *  (w.,w.)  f  H  (D)  x  H  (D)  and  the  linear  functional 


F(W)  =  (f  w  +  f  w  )  ds* 

J3fl  1 

The  weak  solution  U  €  H^D)  x  (D)  satisfies 

(7.4)  B(U,W)  =  F(W)  for  any  W  €  H1  (D)  x  H^D). 

The  strain  energy  functional  G(U)  is 


(7-5)  G(0)  -  ia^w  [B  +  <If>2>+  2v 


3u  3v 
3x  3y 


,  l-2v  ,3u  .  3v.  2,,  , 

+  —  <37  +  }dxdy 


-  B(U,U). 


We  will  consider  problem  1  and  problem  2  when  the  imposed  tractions 
lead  to  the  symmetric  and  antisymmetric  mode  of  the  stress  intensity 
solution.  For  both  modes  i  *  1,2  the  solution  has  singular  behavior  at 
the  tip  of  the  crack  (see  [19]). 

1/2 

u±  %  r  ^  $  ±(9) 

Vi  ~  rl/2  ^l,i(9) 

3/2-r  9  9  1 

and  hence  u  ,v  €  H  (Q)  ,  e  >  0,  and  u  ,  v  €  ~  <  S  <  1. 

I  1  p  2 
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The  solution  of  Problem  1  has  to  be  augmented  by  the  conditions 


1^(0, 0) 


v  (0,0)  -  0 


fixing  the  rigid  body  motion.  For  the  solution  of  Problem  2  we  impose  the 


conditions 


u 2(0,C)  =  v2(0,0)  -  v2 (-1,0)  =  0 


fixing  the  rigid  body  motion.  The  strain  energy  of  true  solutions  for 

E  =  1.0  and  v  =  0.3  is 


G(U  )  -  0.6017796916, 


G(U2)  =  0.2370646876. 


By  UQ,  i  =  1,2  we  denote  the  exact  solutions  and  by  Up^,  i  =  1,2, 
the  finite  element  of  Problem  1  and  Problem  2.  The  error  of  the  finite  element 
solution  will  be 


(7.6) 


ui  -  U* 

0  FE 


The  energy  norm  of  the  error  HeJI  is  directly  related  to  the  strain  energy 
of  the  exact  and  finite  element  solution 


(7.7) 


'ei»E  •  G(US)  -  GWFE)  -  -  B<Ure'°FE)) 


and  the  relative  error  in  energy  norm  is  defined  as 


(7.  8) 


"illE,R 


G(Uq)1/2 


We  used  geometric  meshes  with  the  factor  a  for  the  p  and  h-p  version 
and  studied  the  performance  of  two  types  of  meshes  A  and  C  shown  in 
Fig.  7.2  and  Fig.  7.3. 


(T*  0^0^  O 

MESH  C5 


Figure  7.3 


MESH  C6 


Mesh  C 


Let  us  mention  that  the  size  H  .  of  the  smallest  element  is  C  . 

mm 

-4 

For  a  =  .15  we  get  h  .  =  .759  x  10  and  for  a  =  .08  we  get  h  .  - 

mm  s  min 

.327  x  io  5  with  h  =1.  The  ratio  between  the  size  of  the  largest  and 

max 

smallest  element  is  hence  more  than  10^.  The  computation  has  been  performed 
by  the  PROBE  program  with  uniform  degree  of  elements  in  double  precision  on 
Appolo  420  (work  length  of  double  precision  is  15  decimals) .  The  computation 
for  the  h-version  was  implemented  by  the  program  FEARS  on  UNIVAC  1100  (in 
single  precision) .  By  N  we  denote  the  number  of  degrees  of  freedom. 


As  indicated  in  the  previous  chapter  there  are  some  constant  b  and  C 
independent  of  N  such  that  for  the  h-p  versions  we  have 


(7.  9) 


6i  'E,R  ~ 


C. 

_ l 


-b^N 


1/3 


i  =  1,2. 


0  E 


Table  7.1  and  Table  7.2  show  the  relationship  between  ||e  ||  ,  and 

1  E 

N,  p,  n,  on  mesh  A  with  a  =  .15  and  y  =  1  (i.e.,  p  =  n)  .  The 

relationships  are  plotted  in  Fig.  7.4  and  7 . 5  on  ln-cubic  root  scale.  The 

curve  of  the  h-p  version  is  almost  a  straight  line  which  is  the  envelope  of 

six  curves  of  the  p-version  for  IS  n  S  6.  This  means  the  asymptotic 

property  is  achieved  already  for  n  *  p  =  2.  The  constants  b^  and  C^, 

i  =  1,2  are  numerically  given  in  the  tables.  For  p  =  6  we  have  b  *  0.670, 
C1  C2 

b„  =  0.668,  and  — i =  1.688,  — ~ —  =  1.306. 

"U0ljE 

In  Fig.  7.6  we  show  the  dependence  of  the  error  on  O.  We  see  that  the 

2 

best  value  of  a  is  close  to  (/2-1)  *  .17  which  is  the  theoretically 

optimal  value  in  one  dimension  (see  [16]). 

Fig.  7.6  shows  the  dependence  of  the  error  on  U  characterizing  the 


relation  between  n  and  p. 


TABLE  7.2 

Relationship  between  jje2||E  R'N'n'P'b2  and  C2  for  ***  h-p  Version 
Problem  2  (E=l,v=0.3)  on  Mesh  A,c=0.15,y=l 


Mesh 

B 

N 

„l/3 

||e2 Je/ R  % 

b2 

rTTT~ 
3  IuoIIe 

*1 

B 

8 

2 

43.74 

.626 

1.664 

*2 

2 

47 

3.61 

15.97 

.742 

1.781 

*3 

3 

120 

4.93 

5.91 

.772 

1.592 

E9 

B 

255 

6.34 

2.02 

.718 

1.395 

*5 

5 

476 

7.81 

.71 

.668 

1.309 

A6 

6 

807 

9.31 

.26 

.668 

1.306 

100  200  400  600  800  1000  1200  1500 


Figure  7.6 

Relative  Error  in  Energy  Norm  vs.  Degrees  Of  Freedom  for  the  h-p  \fersion 
Problem  2  (E=l,v«.3)  on  Mesh  A  ,l«;n<6 


Figure  7 . 7 

Relative  Error  in  Energy  Ilona  vs.  Degrees  of  Freedom 
for  the  h-p  version 

Problem  2  (E*l,  v=.3)  on  Mesh  A  ,  l<n<6,  o».15 

n 

1  p=n- 1,  u**8~1.0,  2  p-n,  u*l, 

A  p«n+l,  y-1.0-1.2,  3  p-n+2 ,  p-1.2,-1. A, 

5  p"n+  3,  v»l.A“1.6 
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Table  7.3  gives  the  values  of  and  C^.  We  see  that  the  value 

p  =  1  is  the  best  one  and  the  sensitivity  of  the  error  on  p  is  not  large. 
The  results  for  Problem  1  and  Problem  2  are  very  similar.  Problem  1  has 
been  computed  with  meshes  A  and  C.  To  achieve  the  same  accuracy,  mesh  C 
required  more  computer  time  than  mesh  A. 

It  has  been  observed  and  rigorously  analyzed  in  [8], [22]  that  the 
p-version  is  insensitive  to  Poisson  ratio  v.  Numerical  results  show  that 
the  h-p  version  is  insensitive  to  change  of  Poisson  ratio  v.  In  Fig.  7.8 
the  curve  for  v  =  0.49  is  almost  parallel  to  that  for  v  =»  0.3* 


The  comparison  of  any  codes  is  a  very  delicate  question  because  of  the 
aims  of  computation,  the  reliability  of  the  computed  results,  the  ratio  between 
the  human  and  computer  cost  in  the  project,  etc. 

We  address  here  crudely  the  cost  of  achieving  the  same  accuracy  by  the 
h-version  (elements  of  degree  1)  with  nearly  optimal  mesh  and  h-p  version 
with  nearly  optimal  mesh  and  degrees  of  elements.  The  comparison  is  based  on 
the  computation  of  Problem  1  program  FEARS  and  PROBE  (see  Fig.  7.9). 

The  error  of  the  h-version  with  the  optimal  mesh  is  asymptotically  (p  =  1) 


|e||E  *  C2N 


and  for  the  h-p  version  we  have 


tv  C^e 


The  computer  cost  for  the  h-version  is  roughly 


W  -  C  +  DN 


1 

0. 8-1.0 

1.0-1 . 2 
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n 
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with 


1.5 


FEARS  uses  elimination  method  and  a. 


1.9.  For  the 


^  «»**»«  v>ww  VJ.ilUJ-‘*k*U4.WU  UlVrbUWU  ttUU  ^ 

h-p  version  in  our  case  the  total  number  of  arithmetic  operations  is 
asympto  tically 

7/3 

W  =  C  +  DN  ' 


1/3  2  2 /  3 

when  N  independent  matrices  (condensation)  of  order  p  =  0(N  '  )  are 

decomposed  and  additionally  a  band  matrix  of  size  N*"'  and  band  0(N  /J) 

is  decomposed.  PROBE  uses  the  front  solver.  The  cost  of  computation  of  the 

7/3 

microstiffness  matrices  is  also  0(N  ).  The  main  cost  is  in  the  data 

management . 

Assume  that  effectively 

ah 

W  ~  N 


for  the  h-version  with  optimal  mesh  and 


W 


N 


h,p 


for  the  h-p  version  with  optimal  mesh.  We  can  compute 


In  N  (e) 

*<e)  '  m  h.  U) 

h,p 


which  is  the  ratio  between  and  3^  leading  to  the  same  accuracy  for 

given  relative  error  e.  By  N,  (e)  (resp .  N,  (e))  we  denote  the  number 

h  h ,  p 

of  degrees  of  freedom  for  achieving  the  desired  accuracy.  Toble  7.4  shows  the 

accuracy  £,  N,  (e),  Nt  (e)  and  ij;(£)  . 
h  h,p 

The  hierarchic  structure  of  the  elements  leads  to  hierarchic  elimination 
so  that  the  computation  of  lower  degree  elements  (for  fixed  mesh)  is  obtainable 
by  0(N)  operation. 


TABLE  7.4 


Comparison  of  the  Number  of  Degrees  of  Freedom  for  lie  ||  =  £% 

1  E,R 

between  the  h-Version  with  Optimal  Mesh,  and  the  h-p  Version  with 

Mesh  An  .  Un$6,c«.15,u»l  on  Problem  1  (E-l,v-.3) 


e 

30 

20 

10 

5 

3 

1 

m 

Nh(c) 

— 

56 

118 

446 

1397 

3  883 

34  945 

139  782 

Nh.P<£) 

29 

48 

96 

165 

232 

450 

657 

*  <e> 

1.195 

1.232 

1.337 

1.418 

1.517 

1.712 

1.826 

808  1.367 


TABLE  7.5 

Estimated  Error  of  the  h-p  version 
Problem  1  (E=l,v  =  .3)  on  Mesh  , I4n*6 , o=. 15 ,p=1 


ir«nE 


2.9596E-1  2. 9662E-1  60.83  60.92 


2  9.8774E-2  9.8511E-2  20.28  20.23 


3  3.7033E-2  3.7055E-2  7.606  7.611 


1.2489E-2  T  1.2500E-2  2.565  2.567 


5  I  *4.3359E-3  4.3691E-3  *  .891 


E, R  '  Hell  E , R  %  niell£“  HeHE)/l|e|lE% 


.2189 


-.2669 


.0606 


.0926 


.6689 


Table  7.4  shows  that  1%  accuracy  is  very  expensive  to  get  by  the 
h-version  with  the  elements  of  degree  1.  The  accuracy  0.5%  is  probably 
not  achievable  at  all.  The  h-p  version  allows  us  to  use  a  relatively  very 
small  number  of  elements  to  obtain  high  accuracy. 

The  mesh  design,  although  critical,  is  not  too  dependent  on  the 

geometry  provided  one  is  refining  the  mesh  geometrically  with  a  =  .15 

around  every  singular  point. 

In  practice  the  true  solution  Uq  is  unknown,  but  error  measured  in 

the  energy  norm  can  be  estimated  from  the  energy  norm  of  the  finite  element 

2  2 

solution.  Let  =  ||Ugj|E  and  E(n)  =  ||  UFE||  £  on  Mesh  1  2  n5  6. 

By  (7.1)  we  have  for  any  u  >  0  asymptotically 

E(n)  =  Eq  -  Ce~bn. 

Therefore 

A (n)  ■  E(n+1)  -  E(n)  *  Ce  bn(l-e  b)  , 


b 


In 


A  (n) 
A(n+1)’ 


(  -A(n)  )n/(l  -  - ^ 
W+l);n/U  A(n) 


(7.10) 


is  the  a-posteriori  estimation  of  the  accuracy.  Further 


(7.11) 


(E(n)+]|?||2)1/2 


_ 1  _ 1 _ 

a+E£nI)l/2  =  E£nl  ,  _A(n+lKf/2  * 

(  jjelf  (  A(n)  U  A(n)  V 

£ 


Table  7.5  shows  that  the  true  error  in  energy  norm  is  estimated  by 

llellE“ll«llE  . 

Ilellg  with  — ii^jj -  -  1%.  Such  high  reliability  of  the  error  estimation  is 

£ 

achieved  for  every  n  because  the  h-p  version  has  exponential  rate  of 
convergence  for  very  low  n  and  p  and  stably  maintains  this  behavior 
through  all  n. 

Although  the  above  numerical  results  are  in  no  way  exhaustive,  they 
do  suggest  the  following: 

(1)  The  exponential  rate  of  convergence  of  the  h-p  version  agrees  in 
practical  accuracy  range  with  that  predicted  by  the  asymptotic  approximation  theory. 
The  asymptotic  rate  can  be  achieved  by  low  p  and  n,  e.g.  n  =  p  =  1,2. 

(2)  The  coefficient  b  is  related  to  singularity  as  well  as  the 
factor  a  and  y.  It  seems  that  a  ^  0.15  is  an  optimal  one.  Optimal 

y  for  uniform  p  depends  on  a,  $,  d  mesh  as  well  as  N.  The  computation 

on  mesh  A,  ff  ^  .08,  1  2  n  5  6  and  mesh  C  ,  cr  -  .15,  1  2  n  2  6  shows 
n  n 

that  optimal  y  increases  with  N.  If  the  asymptotic  value  y^  of  the 
optimal  y  is  known,  one  should  select  a  value  for  y  less  than  y^. 

3)  As  the  p-version,  the  h-p  version  is  insensitive  to  change  the 
Poisson  ratio  in  plane  elasticity  equation.  For  v  near  0.5,  the 
exponential  rate  of  convergence  can  be  achieved  for  some  low  n  and  p. 

The  curve  of  the  h-p  version  is  shifted  to  the  right  without  affecting  the 
coefficient  b; 

4)  Since  the  exponential  rate  is  achieved  for  low  n  and  p,  the 
a-posterior  error  estimate  coincides  with  the  true  error  very  well  from  low 
p  and  n  to  high  p  and  n. 

5)  If  a  higher  accuracy  is  required  then  only  h-p  version  is  a 


practical  method. 
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•  To  provide  a  limited  consulting  service  in  all  areas  of  numerical 
mathematics  to  the  University  as  a  whole,  and  also  to  government 
agencies  and  industries  in  the  State  of  Maryland  and  the  Washington 
Metropolitan  area. 

•  To  assist  with  the  education  of  numerical  analysts,  especially  at  the 
postdoctoral  level,  in  conjunction  with  the  Interdisciplinary  Applied 
Mathematics  Program  and  the  programs  of  the  Mathematics  and  Computer 
Science  Departments.  This  Includes  active  collaboration  with  govern¬ 
ment  agencies  such  as  the  National  Bureau  of  Standards. 

•  To  be  an  international  center  of  study  and  research  for  foreign 
students  in  numerical  mathematics  who  are  supported  by  foreign  govern¬ 
ments  or  exchange  agencies  (Fulbright,  etc.) 

Further  information  may  be  obtained  from  Professor  I.  BabuSka,  Chairman, 
Laboratory  for  Numerical  Analysis,  Institute  for  Physical  Science  and 
Technology,  University  of  Maryland,  College  Park,  Maryland  20742. 
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